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Abstract. We consider the 2-dimensional focusing mass critical NLS with an 
inhomogeneous nonhnearity: idtu + Au + k{x)\u\'^u = 0. From standard argu- 
ment, there exists a threshold Mk > such that solutions with ||u||£2 < Mk 
are global in time while a finite time blow up singularity formation may occur for 
||m||^2 > Mk- In this paper, we consider the dynamics at threshold || no 11^2 = Mk 
and give a necessary and sufficient condition on k to ensure the existence of criti- 
cal mass finite time blow up elements. Moreover, we give a complete classification 
in the energy class of the minimal finite time blow up elements at a non degen- 
erate point, hence extending the pioneering work by Merle [23| who treated the 
pseudo conformal invariant case k = 1. 



1. Introduction 

We deal in this paper with a two dimensional focusing mass critical nonlinear 
Schrodinger equation with an inhomogeneous nonhnearity: 

fATTQ\ f idtu = -Au-k{x)\u\^u, G [0,r) X 

^^^^^^ \ n(0,x) =no(x), uo:M.^^C, ^^'^^ 

for some smooth bounded inhomogenity /c : — )■ M^. This kind of problem 
arises naturally in nonlinear optics for the propagation of laser beams. From the 
mathematical point of view, it is a canonical model to break the large group of 
symmetries of the k = 1 homogeneous case. 

From standard Cauchy theory and variational argument, [33], there exists a 

critical number > such that for initial data with 

11^011^2 < Mk, 

the unique solution to (jl.ip is global in time and bounded in H^, while on the 
contrary finite time blow up may occur for ||uo||l2 > M^. We address in this paper 
the question of the dynamics at threshold 

11^11^2 = Mk. 

In the pioneering breakthrough work [23j, Merle proves in the homogeneous case 
k = 1 the existence and uniqueness of a finite time blow up solution with critical 
mass. In this case, existence and uniqueness both rely dramatically on the explicit 
pseudo conformal symmetrj{3 which is lost in the inhomogeneous case. Merle initi- 
ated in ^25j the study of the inhomogeneous problem and obtained in particular a 
sufficient condition to ensure the nonexistence of critical blow up elements for the 
inhomogeneous problem 



^see (fL4)l 

^see Theorem ll.ll b elow 
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In this paper, we give a necessary and sufficient condition on the inhomogeneity 
k to ensure the existence of critical blow up elements, and prove the strong rigidity 
theorem of uniqueness of such minimal blow up elements in the energy class despite 
the absence of pseudo conformal symmetry. 

1.1. The homogeneous case. Let us start with recalling some well-known facts 
in the homogeneous case k = 1: 

( idtu = -Au-\u\'^u, G [0,r) X M2, 

\ u{0,x) =uoix), UQ-.R'^^C. ^ ' 

solutions satisfy the conservation of energy, norm and momentum: 
L^norm : J|w(t,x)p = J \uq{x)\'^ , 

Energy : E{u{t, x)) = \ J | Vn(t, x)|2 - ^ / |u(t, x)p+^ = E{uq), 
Momentum : M{u) = Im (J Vwu) = M{uo), 

and a large group of symmetries leaves the flow invariant: if u{t, x) solves (|1.2p . 
then V(Ao, to, xq, /3o, 7o) G M+ x M x x x R, so does 

v{t,x) = xfu{t + to, Xox + xo- Pot)e'^'^''~^^'^e'^°. (1.3) 

A last symmetry is not in the energy space but in the virial space S = {xu G 
L^} nH^, the pseudo conformal transformation: if u(t,x) solves (|1.2p . then so does 

1 fl x\ la^P 

Following [10] ■ [17] ■ let Q be the unique nonzero positive radial solution to 

AQ-Q + Q^ = 0, (1.5) 

then the variational characterization of Q ensures that initial data uq G with 
II^o||l2 < yield global and bounded solutions T = +oo, [33j. On the other 

hand, finite time blow up may occur for data ||iio||/^2 > ||(5||2,2(]g2). At the critical 
mass threshold, the pseudo conformal symmetry ()1.4p applied to the periodic solitary 
wave solution u{t,x) = Q{x)e^^ yields a minimal mass blow-up solution: 

5(t,x) = ^Q(|)e*^-i, ||5(t)||^2 = ||Q||^2 (1.6) 

which blows up at t = 0. In [23], Merle proves the uniqueness of the critical mass 
blow up solution: a solution u G with ||^io||L2 = ||Q||2,2 and blowing up at t = is 
equal to S{t) up to the symmetries of the flow. Through the conformal invariance, 
this results yields a complete dynamical classification of the solitary wave as the 
only nondispersive solution in S with critical mass. 

1.2. The inhomogeneous case. We now come back to (jl.ip . The conservation of 
momentum no longer holds but norm and energy are still conserved: 

L^norm : / |n(t, x)^ = / |uo(x)p. 

Energy: E{u{t,x)) = ^ f \Vu{t,x)\'^ - f k{x)\u{t,x)\^+^ = E{uo), 

N 

The canonical effect of the inhomogenity is to completely destroy the group of 
symmetry ()1.3p R. and in this sense (jl.ip is a toy model to analyze the properties of 
NLS systems in the absence of symmetries. 



'up to phase invariance 
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Prom standard variational techniques and a virial type argument, Merle has derived 
in [25j the criterion of global existence for (jl.ip : given k > 0, let be 

and let 

^2 = max/cfx) < +00, 

then initial data with 

\\uo\\l2 < Mk = WQk^h^ (1-8) 
yield global and bounded solutions while finite time blow up may occur for 
||mo||l2 > Mk- Moreover, the localization of the concentration point and a sufficient 
condition for the non existence of critical blow up elements are obtained: 

Theorem 1.1 (Merle, [25j). Assume that k belongs to C^(M^). 

1. Localization of the concentration point; Let u £ with \\u\\i^2 = Mk he a 
solution to p.ip blowing up at T = 0. Assume that {x G M? such that k{x) = 

is finite. Assume also that there is 6 > and R > such that k{x) < k2 — S for 
\x\ > R. Then there is xq G such that k{xo) = ^2 o-nd 

\u{t)\'^ \\Qk{xo)\\h^x=xo as t-^0. 

2. Criterion of non existence; Assume that k{xo) = ^2 an-d 

\7k{x) ■ (x — xq) < —\x — xqI^^'^^ near xq and for some > 0, (1-9) 
then there are no critical mass blow up solution at xq. 

In other words, the concentration of a possible minimal blow up solution must 
occur at a point xq where k reaches its maximum, and hence in particular: 

\/k{xo) = 0, 

and the repulsivity condition (|1.9p in fact implies the non existence of such an ob- 
ject. Note however that (|1.9p implies that k is not at xq, and hence Theorem ll.il 
leaves completely open the case of a smooth k which we address in this paper. 

In the presence of a very smooth and flat k at xq, the existence of critical elements 
may be derived using a brutal fixed point argument, [25], |24) . [6], see also |14j for 
related statements in the setting of nonlocal nonlinearities. This argument has 
been recently sharpened by Banica, Carles, Duyckaerts |2j for (jl.ip by adapting 
the method designed by Bourgain and Wang [J] and further revisited by Krieger 
and Schlag [16]: after linearizing the problem close to the explicit S{t) approximate 
solution, one uses modulation theory and energy estimates to treat perturbatively 
the unstable modes and integrate the system backwards from the singularity 0. This 
allows one to lower the flatness of k and indeed the existence of critical mass blow 
up solution is obtained under the flatness assumption: 

Vk{xo) = V^k{xo) = 0. 

Let us say that in this approach, the problem is treated perturbatively from the 
homogeneous case and xq can be taken to be any point where k is flat enough. The 
case of smooth k with nondegenerate Hessian at xq i.e. 

V^/c(xo) < 



'similarly like in scattering theory 
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is out of reach with these techniques because it generates large deformations which 
live at the same scaling like S{t) as we shah see later. 

Let us also say that even when existence is known, uniqueness in the energy class 
as proved in [23j for the homogeneous case does not follow. See [1] for a pioneering 
investigation of this delicate problem in the setting of a domain. Here the strength 
and the weakness of Merle's pioneering proof for = 1 is that it fundamentally 
relies on the pseudoconformal symmetry ()1.4p which is lost here. The challenge is 
hence to provide a more dynamical proof of the classification theorem in the absence 
of symmetry. Eventually, let us recall that uniqueness is one of the keys towards the 
derivation of chaotic continuation after blow up results, [25], which is a canonical 
open problem for inhomogeneous problems like (jl.ip . 

1.3. Statement of the result. Let us now fix our assumptions on k where we 
normalize without loss of generality the supremum k2 = 1- 

Assumption (H): k is D W^'°° with 

< ki < k(x) < 1 and maxk(x) = 1 is attained. 

From (|1.7p . (jl.Sp . the critical mass is then 

Mk = \\Qh2. 

We start with the derivation of necessary conditions for the existence of a critical 
mass blow up solution which is based on a slight refinement of the variational 
techniques introduced in |25| : 

Proposition 1.2 (Necessary condition for the existence of a critical blow up el- 
ement). Let u with ||m||^2 = ||(5||l2 be a solution to (jl.ip which blows up at time 
T = 0, then there exists xq G M? such that 

k{xo) = 1 

and u blows up at xq in the sense: 

kWI' ^ IIQIIi24=xo «s i^o. (1.10) 

Moreover, the energy Eq of u satisfies: 

Eo + ^lv\{xo){y,y)Q^>0. (1.11) 

In order to classify critical mass blow up solutions, we hence pick xq G such 
that k^xo) = 1 and focus in the whole paper onto the case of a nondegenerate 
Hessian 

V^k{xo) < 

which we expect to be the most delicate one. We claim that the lower bound (jl.lip 
is sharp, and the following theorem is the main result of this paper: 

Theorem 1.3 (Existence and uniqueness of a critical element at a nondegenerate 
critical point). Let xq G with 

k{xo) = 1 and V^k{xo) < 0. 

Then for all Eq satisfying (jl.lip . there exists a unique up to phase shift critical 
mass blow up solution to (|l.ip which blows up at time T = and at the point xq in 
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the sense of (jl.lOp . with energy Eq. Moreover, 



lim/m VnuJ = 0. (1.12) 

In other words, for k smooth, there exists a critical mass finite time blow up so- 
lution if and only if the supremum of k is attained, and the corresponding minimal 
blow up elements at a non degenerate blow up point are completely classified. 

Comments on Theorem 11.31 

1. Structure of the minimal blow up elements: Let 

Eo = Eo + l[v'kmy,y)Q'>0, Co = (1.13) 

and 

Sc,(t,x) = (^^^^) (1.14) 



be the exact pseudo conformal blow up solution to the homogeneous problem (|1.2p 
with modified energy p.l3p . then up to a fixed phase shift, the critical mass blow 
up solution of Theorem 11.31 decomposes into a regular and a singular part 

u{t) = Scoit) + u{t) 

with 

u{t) ^0 in H'^ as t ^ 0. 
In particular, the blow up speed is given by the conformal law: 

|Vn(t)U2 ~ (1.15) 

There are two fundamental remarks. First, the manifold of critical mass blow up 
solutions at a nondegenerate point is now strictly smaller than in the homogeneous 
case and is submitted to the constraint of the non trivial lower bound on the energy 
(jl.lip and an asymptotic zero momentum condition ()1.12p . This non trivial struc- 
ture is very much due to fine non trivial effects induced by the inhomogeneity k. In 
view of ()1.14p , the leading order effect of k is to induce a pseudo conformal type blow 
up but parametrized by the modified energy (|1.13p . This shows in particular that 
the case V^/c(xo) 7^ cannot be treated perturbatively from the homogeneous case 
as for a given initial data -and hence a given level of energy-, the blow up profiles of 
minimal mass blow up solutions for (jl.ip and ()1.2p are strictly different due to the 
energy shift (|1.13p . Second and most importantly, the fact that the blow up speed is 
still given by the pseudo conformal law (|1.15p is very much a surprise. Indeed, the 
nonvanishing condition V^A:(xo) 7^ induces deformations which leave at the same 
scaling like the leading terms driving the homogeneous pseudo conformal blow up. 
We shall see that the derivation of the law (jl.lSp relies on an algebraic cancellation. 
Both these issues are detailed in the strategy of the proof below. 

2. The degenerate case: We expect the methods presented in this paper to also 
be able to treat the degenerate case Vk^xo) = \7'^k{xo) = which should in fact be 
easier to handle. In particular, this becomes a perturbative problem with respect 
to the homogeneous case and we expect the manifold of minimal mass blow up 
solutions to be of the same size like in the homogeneous case. More precisely, we 
conjecture the following: for all Eq > and Mq G M?, there exists a unique up to 
phase shift critical mass blow up solution to (jl.ip which blows up at time T = 
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and at the point xq in the sense of (jl.lOp . with energy Eq and asymptotic kinetic 
momentum 



Note that the existence part for Mq = is proved in [2]. 

3. More blow up solutions: The existence of stable log-log type blow up dynamics 
for super critical mass initial data HuoHl^ > IIQIIl^ can be derived by adapting the 
argument in |29j, see | i30j . ^ for related results. Such solutions would blow up with 
speed 



This is a consequence of the fact that in the log-log regime, the inhomogeneity k can 
very much be treated perturbatively and the problem be essentially reduced to the 
homogeneous one. This is a manifestation of the stability of the log-log blow up, 
|31j . Theorem 11.31 shows that the structure of critical mass blow up solutions is of 
course much more fragile due in fact to its intrinsic instability. In the homogeneous 
case, there also exist excited pseudo conformal blow up solutions corresponding to 
the infinite sequence of solutions to (jl.Sp in dimension N > 2, which display 
a nondispersive behavior in the sense that all the mass is put into the singularity 
formation: 



The existence and possible uniqueness properties of such objects with higher 
norm is open, see [2g'j for some classification results about such solutions in the 
homogeneous case. 

4- Extension to higher dimensions: We have decided to restrict the analysis to 
dimension N = 2 which is the physically relevant case. The case = 1 could 
be treated completely similarly. For higher dimensions, the only restriction of our 
analysis is of technical nature and one would have to face the problem of the lack 
of differentiability of the critical nonlinearity ii|n|^ for large A^. 

5. Connection to previous works: This work lies within the question of classifica- 
tion of critical elements. This kind of question is related to Liouville type theorems, 
|12j , and is connected to the question of the dynamical classification of solitary wave 
solutions which is only at its beginning, see [ 1 9j , |13) . In particular, our techniques 
borrow from [19\ the energy estimates type of strategy, and are in spirit related 
to the recent work by Duyckaerts and Merle |Sj where deep rigidity properties of 
critical dynamics of the energy critical homogeneous NLS are proved. One main dif- 
ference with these works however is the fact that we are dealing here with a slowly 
polynomially converging problem, while |8], |19) deal with exponential decay which 
somehow allows a more clear decoupling of the interactions. 

Eventually, let us say that our strategy both for existence and uniqueness does 
not rely at all on the existence of the pseudo conformal symmetry for the homo- 
geneous problem. On the contrary, the singular dynamics is extracted directly in 
the renormalized variables and the presence of the symmetry only allows for some 






as t — )• T with u* E L^. 



u{t)\' 




Uo\\l'2^x=xo- 
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-unnecessary- algebraic simplifications. Let us recall that the striking example of 
the L?' critical KdV shows that critical mass blow up dynamics do not always exist, 
[21| . and non uniqu^ critical mass blow up dynamics are exhibited in \Z\ for the 
critical Keller-Segel model. We expect our strategy to be robust enough to provide a 
canonical framework to address the question of existence and uniqueness of minimal 
blow up dynamics for various nonlinear dispersive equations in the absence of any 
explicit pseudo conformal symmetry. 

Acknowledgments: Both authors are supported by the French Agence Na- 
tionale de la Recherche, ANR jeunes chercheurs SWAP. P.R is also supported by 
ANR OndcNonLin. 

1.4. Strategy of the proof. Let us give a brief insight into the strategy of the 
proof of Theorem 11.31 



step 1 Refined blow up profile and formal derivation of the laws of motion. 

Let us start with assuming the existence of a critical mass blow up solution at 
T = and xq = and obtaining dynamical informations on such a solution. From 
standard variational argument, u must admit near blow up time a decomposition 
of the form: 

x) = / {Q + e) (t, ^ ~ "^-^^ ) e^'^(*^ with \\e{t)\\Hi ^ as t ^ 0, 
X-{t) V Ki) J 

see [33], [25]. In fact, following [26], [27], [31], this decomposition should be sharp- 
ened by deforming the main nonlinear part in the direction of pseudo conformal 
blow up. More precisely, in the homogeneous case, the explicit pseudo conformal 
blow up solution (jl.6p should be thought of as being the following object: 



u{t,x) = -i— Q,(,)^^(,) (t, ^^^) e^^W with Q,Ay) = Q{v)e-"-^e'P-y 
\~{t) \ ^V'l J 

and where the geometrical parameters (A(t), are driven by the dy- 

namical system: 

bs = -b\ -y = ^ y = 2/3, /3s = -6/3, 7s = -1 (1.16) 
where s is the global rescaled time 

ds__ ^_ 

Integrating the system with boundary condition A(0) = -corresponding to blow 
up at t = 0- and a(0) = -corresponding to blow-up at xq = 0- involves three 
integration constants which can de directly related to the energy and momentum of 
the initial data, and the choice of a fixed phase shift. 

The first step in the proof of Theorem 11.31 is to similarly decompose the solution as 
u{t,x) = -^{Qp + e) (t, e-^W with P = (6,A,/3,a) (1.17) 

where the profile Q-p takes fully into account the deformation induced by the k 
inhomogeneity, and is an approximation of high order -0(A^)- of the renormalized 
equations. This construction generalizes the one introduced in |27) . see also [15j for 



infinite time blow up 
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a related approach in a different setting. Under the assumption that the excess of 
mass £ is of lower order in some suitable sense, one may formally extract the leading 
order ODE driving the geometrical parameters which roughly takes the form: 

bs = -b^ /3, = -6/3 + ci(a)A, y = y = 2/3, K*) ~ ^ (1-18) 

where ci{a) is a linear form in a. The last two equations are standard. The 
first equation for b is surprising as inspection reveals that we should expect from 
V^kiO) / 0: 

bs = -b^ + coX\ 

and A ~ 6 in the conformal regime. It is a specific algebraic cancellatioi^ which 
leads to cq = 0. The second law for /3 reflects the difference between the degenerate 
case for which ci = and the non degenerate case for which ci(a) ~ a. In this last 
case, the reintegration of the scaling law yields 



X{s ^ 



s ' ' s 



On the other hand, from (|1.17p . 

jr^Im(^j Vnn^ ~ Mq (1.19) 

and hence if the asymptotic momentum Mq is non zero asymptotically, the reinte- 
gration of the translation parameter yields 

at ~ Mq and thus a(t) ~ Mo|t| ~ — -. 

s 

Injecting this into the momentum law of /3 yields: 

=i[/3. + 6/3]~ci(a)~^. 

This generates a logarithmic divergence in s by integration which contradicts (|1.19p 
for Mq 7^ 0. Such a phenomenon is absent is the degenerate case V^A;(xo) = for 
which ci = 0, and hence the asymptotic value Mq becomes a free parameter when 
integrating (jl.lSp as in the homogeneous case. 

step 2 Construction of a critical mass blow up solution. 

The generalization to a critical setting of the subcritical strategy of slow variable 
ansatz developed in |15| allows us to derive an approximate critical mass solution to 
a sufficiently high order. The derivation of an exact minimal mass blow up element 
then follows from a standard compactness argument using energy estimates and 
integrating the flow backwards from blow up time as initiated in [23] . see also |19) . 
[22], [15| . Indeed, one decomposes the flow as 

U{t, X) = Uapp{t, x) + U{t, x) 

where u roughly solves: 

idtu + Lxu = l.o.t. (1.20) 

Here the linearized Hamiltonian Lx is after renormalization the matrix linearized 
operator close to the ground state L = (-L_|_, L_) with: 

L+ = -A + l-3Q^ L_ = -A + 1-Q2. (1.21) 



6, 
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Similarly like in [32], a direct energy estimate on (|1.20p however fails and generates 
quadratic errors due to the time dependence of the Hamiltonian which cannot be 
treated perturbatively. We hence modify the energy functional by adding a local 
Morawetz type term and derive a monotonicity formula of the form: 

|2 



^|l|Vn||i2 + ^^ + 9(^^ 9,nn^| > + O(Lo.t.), (1.22) 

which for a sufficiently small error in the constructed approximation is enough to 
close the existence part. 

step 3 Derivation of the blow up speed assuming dispersion. 

We now turn to the proof of uniqueness. Standard variational arguments ensure 
that a critical mass blow up solution must admit a decomposition of the form (jl.lTp , 
but there holds no a priori upper bound on the blow up rate in this regime. In fact, 
if we rewrite ()1.17p in the form 

„(,,)._J_«,(,i_^V-"" + ^. 

then standard variational estimates yield: 

limsup \\u{t)\\Hi < C, lim \\u{t)\\L2 = 0. 

The key to the derivation of the sharp blow up speed is to improve this energy 
bound for a dispersive estimate: 

u{t) ^ in H^. (1.23) 

The proof of ()1.23p is the heart of the paper and relies on the derivation of a local 
Morawetz monotonicity formula which generalizes in some sense (jl.22| ^ and provides 
another Lyapounov type rigidity in the problem in the form 

V£tp+ / luM • (1-24) 

Note that remarkably enough, the quadratic form involved in the above coercivity 
estimate is the linearized quadratic form of the energy for which the coercivity is 
known from the variational characterization of Q, and no further coercivity property 
will be needed in the papei@. 

The control of the local interactions with the ground state part of the solution 
provided by (|1.24p coupled with the critical mass conservation laws will imply all 
together the dispersion (jl.23p and the derivation of the blow up speed: 

6 ~ A ie A(t) ~ \t\. 

step 4 Uniqueness. 

Once the blow up speed is known, one can use the energy machinery p.22p again 
to compare the flow with the constructed solution Uc and derive a bound 

\\u-u,\\hi <\t\'^ (1.25) 

for some constant A. Uniqueness then follow by integrating the flow close to Uc 
backwards from infinity using a dynamical system type of argument and the a priori 




which is useless until we know the blow up speed 
^as opposed to the works [20) . |26) . [S] 
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bound (jl.25p . The difficulty is to beat tlie polynomial growth of the propagator e**^ 
generated by the null space of the -non self adjoint- operator L in a context where all 
possible gains are also only polynomial^. Recall indeed that the group of symmetries 
of the homogeneous case generates polynomially growing solutions to the linearized 
flow, explicitly: 

L+(Vg) = 0, L^{yQ) = -2VQ, (1.26) 

L_Q = 0, L+(AQ) = -2Q, L_(|y|2g) = -4AQ, L+p=\y\'^Q, (1.27) 

where p is the unique radial solution to L+p = \y\'^Q- Beating the induced 
a priori growth of the propagator e**^ requires a sufficiently large a priori constant 
A = 3 in (|1.25p which in turn forces the construction of an approximate solution to 
a sufficiently high order. 

The rest of the paper is as follows. In section O we construct an approximate 
solution to the order 5 in renormalized variables, Proposition 12.11 In section [3l 
we derive a mass critical energy identity. Lemma 13.31 which is enough to prove 
the existence of minimal mass blow-up solutions. Proposition 13.61 In section 31 we 
start the proof of uniqueness. We derive the local virial identity. Proposition 14. 3| 
which is the key to the control of the blow up speed and dispersion, Proposition 14.61 
In section [5l we iterate the use of the energy identity of Lemma 13.31 to control the 
growth induced by the null space. Lemma [5. 31 and conclude the proof of uniqueness. 

Notations We let 

if,9) = / f{x)g{x)dx 
be the scalar product. We introduce the differential operator 

A/ = / + y • V/ 

and recall that 

(/,Ag) = -(A/,g). 



2. Approximate solutions and leading order dynamics 

From now on and for the rest of this paper, we assume that k satisfies As- 
sumption (H). Moreover, without loss of generality, we assume that k attains its 
maximum at xq = which is nondegenerate: 

A;(0) = 1, Vfc(O) = 0, V^fc(O) < 0. (2.1) 

Our aim in this section is to construct a high order approximate critical mass blow 
up solution by following the slow variables method implemented in particular in 
|15| . The outcome is the derivation of the leading order dynamical system (|1.18p 
driving the geometrical parameters attached to the refined modulation theory. In 
particular, we will see how an algebraic cancellation leads to the preservation of the 
pseudoconformal speed, and a new rigidity occurs for the motion of the center of 
mass. 



'and not exponential as in [19| . [S] where a similar issue arises 
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2.1. The slow modulated ansatz. Let us consider the general modulated ansatz: 



which maps the finite time blow up problem (jl.ip onto the global in time renormal- 
ized equation (|2.2p : 

k{\{t)y + ait)) , ,2 
k[a{t)) 

• As A -cts XVk(a(t)) \ ^ 

with 7s = 7s — 1. We freeze the homogeneous laws 

Aq . a. 



and let 



so that (12.21) becomes: 



^ ^, ^ 2/3, 7. = l/3|^ 



w{s,y) = v{s,y)e^^~^^-y, (2.3) 



z^su; + A..; - + (6, + 62)k5^! _ |(^^ + 6^) . y + iA/3 • ^^^ll^l Tx; 

K(a(t)) 



(2.4) 



Note that Q itself provides an approximate solution of order but this is not enough 
to carry out the analysis and to exhibit fine effects induced by the homogeneity k. 
In order to construct a higher order approximation, we adapt the strategy designed 
in [15) in a subcritical setting to the critical case. Given small parameters 



P = (6,A,/3,a), 
we look for a solution to (12.41) of the form 



w{s,y) = Pj,(^,){y), y = -&, y = 2/3, Is = bs + b^ = 0, ps+b(3 = 13{X,a). 

This maps (|2.4p onto: 

-ib'^dbPv - iXbdxPv + 2ipXdaPv + i{-b/3 + B)dpP'p 

■ li^} ------ 

The profile P-p is computed from an asymptotic expansion near Q and at each step 
the laws for the modulation parameters B which contain the deformation induced by 
the inhomogeneity k are adjusted to ensure the solvability of the obtained equations. 
Indeed, recall that the kernel of the linearized operator close to Q is explicit: 

Ker{L+} = span{VQ}, Ker{L^} = span{Q}, (2.5) 

see |33) , [7] . The 2 instability directions generated by the kernel of are precisely 
adjusted by computing the two parameters B GM?. The instability direction gener- 
ated by the kernel of L- should in principle be adjusted by modulating on the law 
of b. However, a spectacular cancelation (see (j2.16p ) will take care of the kernel of 
L_, which in turn will allow us to keep the same law for b as in the pseudocoiiformal 
regime, i.e. bg + b'^ = 0. 
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Proposition 2.1 (Approximate solution). Let Cq > be a given constant. There 
exists a universal constant c > and a small constant 7/*(Co) > such that for all 
\V = {b, A, /?, a)| the following holds true. There exist B of the form: 

B = Aco(a) + /33A3 + /34A^ (2.6) 

where cq is a linear map on and /33, /34 are vectors in M? , and smooth well localized 
profiles (Tj, Sj)i<j<4 homogeneous of degree j inV, such that 

4 

Pv = Q + Y.^Tj + iSj) (2.7) 
i=i 

is an approximate solution to ()2.5p in the sense: 

-ib^diP-p - iXbdxPv + 2if3XdaPv + i{-b^ + l3)d^Pr 

B.y + ^Xp.^^]pr + APr-Pr 
k{X{t)y + a{t)) p IP |2 . ...^ 

+ M«(t)) ^-'^^^^^ ^'-'^ 

^•p(y) polynomial in V and smooth and well localized in y: Va € N^, 

d^^viy)] < (v''+V{a^ + P^)+(^b--^^V^^e'^-\y\. (2.9) 

Remark 2.2. Remark that the remainder (|2.9p is not of order five with respect to 
aU parameters. This simpHfies the construction and is not a problem since we will 
obtain better decay properties for a, (3 and b — X/Cq. In particular, we will show 
a = 0{V^), (3 = OiV"^) and b- X/Cq = 0{V^) where Cq will be explicitly chosen 
according to (|1.13|l ^1. 

Proof of Proposition 12.11 

The proof proceeds by injecting the expansion (|2.7p in (|2.8p . identifying the terms 
with the same homogeneity, and inverting the corresponding operator. Let us recall 
that if L = (L+,L_) is the matrix linearized operator close to Q given by (jl.2ip . 
then the explicit description of the kernel (|2.5p ensures the standard uniform elliptic 
estimates: 

V/ E (VQ)^, ||e^lJ^lL;Vk2 < ||e2^l^l/||^2, (2.10) 

V/ G (Q)^, ||e*lL:V||H2 < ||e2^l^l/||i2. (2.11) 
step 1 Derivation of the law for the center of mass. 

We inject (j2.7p in (j2.8p and sort the terms of same homogeneity. 

order 1: From VA;(0) = 0, we obtain 

L+{Ti) = 0, L^{Si) = and hence Ti = 0, Si = 0. 
order 2: We identify the terms homogeneous of order 2 in (|2.5p and get: 



L+(r2) = V^k{0){a, y)XQ^ + ^V^k{0){y, y)Q^ - Aco(a) • yQ, (2.12) 

L_(52) =0. 



-'^'^see Remark 
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For S2, we may choose: 

In view of the first equation in (|2.12p and (j2.10p . we may solve for T2 if we adjust: 

V^m{a,y)XQ^ + ^V^k{0){y,y)Q^ - Aco(a) • yQ,VQ^ = 0, 
which computes the value of cq: 

HQ' ) 
HIQ 

step 2 Derivation of pseudo-conformal speed: an algebraic cancellation. 



co(«)), := fl^V''kiO)ie,,a) for j = 1,2. (2.13) 



Recall that in the pseudoconformal regime A(t) ~ t, there holds 

6 = -y = -AAt ~A. 
Equivalently, the law for b is in this case 

bs = -b\ 

Hence the term ^V^A;(0)(y, y)Q^ in the RHS of the T2 equation (j2.12p may poten- 
tially induce a drastic modification of the b law of the form: 

bs = -6^ + CA2 + . . . 

This explains why the non degenerate case is a much larger deformation than 
the degenerate case V^A;(0) = 0. We now claim that in fact the choice C = is 
dictated by the algebraic cancellation ()2.16p . To see this, we compute the system 
at the next order: 

order 3: We identify the terms homogeneous of order 3 in (|2.5p . We get rid of 
the terms Vo? ^ Va/3 and which are remainder terms. We obtain: 

L+(T3)= v''k{0){y,y,y)^Q^ + V^k{0){y,y,a)^Q' - l3sX^ -yQ, (2.14) 
L_(S3) = -XbdxT2 + 2l3XdaT2. 
In view of the first equation in (|2.14p , we may solve for T3 if and only if: 

V^miy, y, y)^Q^ + V'k{0)iy, y, a)^Q^ - P^X' • yQ, VQ j = 0, 
which is true provided we choose /S^ to be: 

W3h = ^J^5^ for J = 1,2. 

In view of the second equation in ()2.14p . we may solve for 5*3 if: 

{-XbdxT2 + 2(3XdaT2, Q) = 0. (2.15) 
We now integrate by parts using Lj^{KQ) = —2Q and (I2.12p to compute: 

- 2(r2,Q) = {L+T2AQ) = ^{V^m{y,y)Q\KQ) = 0, 

which shows that (|2.15p does actually hold. Here we used the spectacular algebraic 
cancellation: 

(y,■y^Q^AQ) =0 fori,/ = 1,2. (2.16) 
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This cancellation is the reason why we did not need to adjust the law for b at the 
level and hence why the critical blow up will still be of conformal type. 

step 3 Higher order corrections. 

At this stage, we have computed the leading order terms in the derivation of 
the modified modulation equations. We however keep expanding the profile to get 
higher order corrections, that will be required for the final step of the proof of 
uniqueness, see section [5l 

order 4: We identify the terms homogeneous of order 4 in ()2.5p . We get rid of the 
terms V^a and P which are remainder terms. Since terms of type (6 — \/Cq)V^ 
are also remainder terms, we may replace everywhere b by A/Cq. We obtain: 

= -/34A^ • 2/Q + /4A^ 

- 

- --pr(J\J^'i, 

where /4 depends on T2, Ts, S'^ and is a function from to M. In view of the first 
equation in (I2.17p . we may solve for T^^ if and only if we have: 

/?4-2/Q + /4,VQ^ =0, 

which is true provided we choose /34 to be: 

_ 2(JMQ) 

(P4jj (7q2^ J - i> ^• 

In view of the second equation in (|2.17p . we may solve for ^4 if and only if we have: 

-^9Ar3,Q^=0. (2.18) 

Since Lj^{KQ) = —2(5, this is equivalent to: 

A2 




2^^(9AL+(r3),AQ) = 0. (2.19) 

Now, in view of the first equation of (|2.14p and the cancellation (|2.16p . we have: 

(L+(r3),AQ) = 0. (2.20) 
In view of (I2.19P and (j2.20p . (j2.18p holds. This concludes the construction of 84^. 



step 4 Proof of ([22]). 



The estimate ()2.9p now follows from the above construction, the uniform exponen- 
tial bounds ()2.10p , (|2.1ip and the explicit polynomial development of the nonlinear 
term ujup. The details are left to the reader. 
This concludes the proof of Proposition 12.11 

2.2. Properties of the approximate profile. Let P-p be the refined pseudo- 
conformal profile given by Proposition 12.11 which has been computed for simplicity 
in the conformal variables (|2.3p . Going back to the self similar variables, we let: 

= Ppe-''^T^+^^-^ (2.21) 
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SO that from direct check: 



Qv 



k{a{t)) 



(2.22) 



where 



(2.23) 



satisfies from (12.9 



(2.24) 



Let us compute the norm and energy of Q-p which wih appear as important 
quantities in the analysis. 

Lemma 2.3 (Invariants of Q-p). There holds: 



J \Qv\' = jQ' + OiV' 



(2.25) 



Proof of Lemma 12.31 

We compute the mass of the approximate solution up to fourth order. Using 
([2:2T]) . we have: 

J \Qv\^ = J{Q + T2 + + 0{r^) = jQ^ + 2(^2, Q) + 2{T3,Q) + 0{V^). 
From (|236]1 : 

(T2,Q) = 0, (2.27) 
while from (|2:20D and L+(AQ) = -2Q, 



(T3,Q)=0, 



(2.28) 
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and ([2:25]) follows. 

For the energy, we have from (|2.2ip : 



E{Q 



r) 



1 f k{Xy + a) 
k{a) 



\Qvf 



(2.29) 



by 



VQ + VT2 + VTa + Vr4 + ( -/3 + ^ ) (^3 + ^4) 



/? - y ) (Q + T2 + T3 + r4) + V53 + V54 



k{Xy + a) 
k{a) 



2\2 



HQ +T2+T3+nr + {83+84)') 



1 1 - \ I Q' + j I \yW + ^-^ I Q' + {T2, -AQ - 



+ (T3,-AQ-Q3) 



( k{Xy 



a 



\ k{a) 



1 Q^ + 0(T'^). 



Injecting the soliton equation -AQ - Q-^ = -Q and 1^^, (l2:28]l into ([229]), and 
finally expanding the inhomogenity k{Xy + a) near using the radiallity of Q yields 

([226]). 

This concludes the proof of Lemma 12.31 



3. Energy estimates and existence of critical elements 

Our aim in this section is to provide a robust framework to prove the existence 
of critical blow up elements in the absence of explicit symmetries. Once the leading 
order behavior is exhibited through the construction of the approximate profile with 
a small enough correction, our strategy is to integrate the flow backwards from the 
singularity. Similar strategies are well known in scattering theory and it was used 
in the (NLS) framework in |24) . [1], |19) . |22) . |15| . The key is energy estimates 
as recalled in the strategy of the proof, with this additional subtlety in the critical 
settin£3 that one needs to add a Morawetz type information to derive a suitable 
Lyapounov type monotonicity information, see ()3.37p . 

3.1. Nonlinear decomposition of the wave and modulation equations. Let 

u{t) G be a solution to (jl.ip on a time interval [^0,^1], < 0. We assume that 
on this time interval, the solution is in what will be proved to be the asymptotic 
regime for critical mass blow up solutions, namely we assume that u{t) admits a 
geometrical decomposition 



with a uniform smallness bound on [tQ,ti]: 

\V{t)\ + Mt)\\m<m<^^- (3.2) 

Moreover, we assume that u{t) has almost critical mass in the sense: Vi G [to,ii], 

|IHIi2-iigiii.|<A^t). (3.3) 



-'^^which was in particular absent in |15| 
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From standard modulation argument, see e.g. [26], the uniqueness of the nonUnear 
decomposition (j3.ip may be ensured by imposing a suitable set of orthogonality 
conditions on e, namely: 

(e2,VS)-(ei,Ve) = 0, (3.4) 
(ei,yS) + (e2,ye) = 0, (3.5) 
-(ei,AG) + (e2,AS) = 0, (3.6) 
(ei,|y|2s) + (e2,|y|'e) = 0, (3.7) 
-(ei,P2) + (e2,/3i) = 0, (3.8) 
where p is the unique even solution to L^p = \y\'^Q and 

Pi + iP2 = p{y)e-'''-^+'^-y. 

Recall the non degeneracy: 

(p,Q) = -\{L+p,KQ) = -l(|y|2Q, AQ) = \\yQ\l2. (3.9) 

These orthogonality conditions correspond exactly in the case V = (0, 0, 0, 0) to the 
null space of the linearized operator close to Q, see ()1.26p . (jl.27p . and the directions 
involved in ()3.4p - ()3.8p provide a first approximation of the null space close to 
From standard argument, the obtained modulation parameters are functions of 
time, see [26j for related statements. Let 

sit) = l ^ (3.10) 

be the rescaled timJ^. then e satisfies for s G [sq, si] the equation: 

(6, + 62)a,s + A (^y + 6^ 5aS (3.11) 

+ (/3, + 5/3 - co(a)A - /SgA^ - (3iX'')dpJ: + ^ (y " ^Z^) 
+ dssi - AMe) + 6Aei - 2/3 (ve, + ^^^^i) " ^'^2 
^ + b^ (AS + Aei) + (7, - |/3|2)(G + £2) 

{bs + b'^)dbe + x(^ + b\dxQ (3.12) 



and 



+ (/3, + 5/3 - co(a)A - /33A^ - /3^X'')dpe + ^ (x " ^^j 
+ dse2 + Mi(e) + 6Ae2 - 2/3 fve2 + ^^^^2) + |/3|'ei 



2 /c(a) 

^ + 5) (AG + Ae2) - (7s - |/3P)(S + ei) 



''^^A slight refinement will be needed, see (|5.105|) . (|5.106p 
-'^'^which will be a global time, s{[to, 0)) — [s{to), +00) 
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where Mi,M2 are small deformations of the linearized operator (L_|-,L_) close to 
Q: 

k{Xy + a) 



Mi{e) 



M2{e) 



-Aei+ei 



-Ae2 + £2 



k{a) 

k{Xy + a 
k{a) 



■((|Qp|2 + 2S2)ei + 2See2), (3.13) 
■((|Qp|2 + 2e2)e2 + 2Seei), (3.14) 



where the nonlinear terms are given by 

Riie) = 3Se? + 2eeie2 + Se^ + |epei, (3.15) 

R2{e) = seel + 2i;eie2 + ^4 + |epe2, (3.16) 

and where ^-p given by (I2.23|) denotes the remainder term in the equation of Qp 
and satisfies (|2.24p . 

Let us collect the standard preliminary estimates on this decomposition which rely 
on the conservation laws and the explicit choice of orthogonality conditions. 

Lemma 3.1 (Preliminary estimates on the decomposition). There holds the bounds 
for s G [so,si]: 
1. Energy bound.' 



62 + |/3|2 + |a|2 + ||e|||,(j^,) [eo + ^J V2fc(0)(y, y)Q^^ +0{V' 
2. Control of the geometrical parameters.- Let the quadratic forms: 



nQW] 



i\y\'Q,p) 



do{a, a) = 7r^7#V^A:(0)(Q, a), di{a, a) = '^''^' do{a, a), 

\\yQ\\L2 4(p,Q) 



(3.17) 



(3.15 



and the vector of modulation equations 
Mod{t) = (bs + b"^ -do{a,a),^s 



|2 + di(a,a),^-2/3,^ 



b, 



/3, + 6/3-co(a)A-/33A3) 
then the modulation equations are to leading order: 

\Mod{t)\ <V^ + V^e\\L2 + \\e\\l, + \\e\\]j, + (^-^^V^ + V{a^ + /J^ 

with the improvement: 



(3.19) 



(3.20) 



+ b 



<V''+VHe\\L2 + 



IL2 



+ Ikll^i + ( ^■ 



+ V{a'^ + (3"^). (3.21) 



Proof of Lemma 13.11 

step 1 Conservation of norm and energy. 

Let us write down the conservation of norm using (j3.ip . (|3.3p : 



j \Q-p + e\'^dy = k{a) j 
from which we obtain: 
23?(e,Q^) + / = k{a 



\u\ 



\u\ 



k{a) / Q'^ + k{a) 



\u[ 



\QT,\^-k{a) / Q 



(3.22) 
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We now use ([2:25]) . (f33|) to conclude: 



2 

Consider now the conservation of energy: 



j Q^^+2ul^j eQv^ + j |ep = 0{V* + V\a\'^). (3.23) 



2 

We expand the nonUnear term: 



(3.24) 



\Qv + e\'^ = |Qp|4 + 4K(e|Qp|2Qp) 

/ 2T? E0 20^ 

+ m{\e\^eQ^) + \e 



We now inject the value of E{Q-p) given by (j2.26p and estimate the cubic and 
higher nonlinear terms using standard Gagliardo-Nirenberg estimates and the a 
priori smallness ()3.2p to derive: 



\Qv?) ' IQ^P V \Qv? 



+ 0(lkllii(R2) + lk||ii(R.)p2) 



+ jft ( £, -AQp - ^^^^^^ + IQpPQp ) + 0(||g|||,i(K^) + T^^^ 



, ^ / 1 - i / ((1 . ^) .f . 4^..., - - ^) ^3 

and hence using (|3.24p : 

step 2 Coercivity of the linearized energy and proof of (|3.17p . 
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We now sum the conservation of mass (j3.23p and the conservation of energy 
(I3:25D . We obtain: 



(3.26) 



/ '^-'^ ((^ + W?) + '1^^^^^ + ^) ^0 
= X^{E, + ^ j V^miy, + 0(||e||3,i + jap + V% 

The remaining Unear term is degenerate from (|2.22p : 

and thus using also the orthogonality conditions (j3.4p and p.6p : 



= 0(p2||e||^,). (3.27) 

lK26ti and (lOT)) imply: 

= A2(^i^o + ^ y"v2A;(0)(y,y)g')+O(lkll?,i+|a|3 + ||e||^ip2^p4). 

We now observe from the proximity of Qp to Q ensured by the a priori smallness 
()3.2p that the quadratic form in the LHS of (|3.28p is a small deformation of the 
linearized energy close to Q. In other words, we rewrite (|3.28p : 

f\ ,2^2 , f-2 V2fe(0)(a,a) 



/l.lV + f jQ^-rm^ /Q^ + l[(L,....) + (^..2,e.)] 

= X^(^Eo + lJ V2fc(0)(y, y)Q^^ + o(||e|||,, + [ap) + 0{\\e\\H.V^ + p4)(3.29) 

We now recall the following coercivity property of the linearized energy which is a 
well known consequence of the variational characterization of Q: 

Lemma 3.2 (Coercivity of the linearized energy, |34) . [26], |29j). There holds for 
some universal constant cq > ; Ve £ , 

(L+ei,ei) + (1^62,82) > co||e|||^i (3.30) 
- i {{si,Qf + (ei, \y\^Qf + {si,yQ? + {s2,P?} • 
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The choice of orthogonaUty conditions (j3.4p - (j3.8p together with the degeneracy 
inherited from (|3.23p : 

|(ei,Q)|2 = o(||e||^i + |a|2) + 0(p4) 

yield: 

(L+ei,ei) + (L_e2,e2) > ^Mm + o{\a\^) + 0{V^). 
Injecting this into ()3.29p and using the nondegeneracy ()2.ip yields ()3.17p . 

step 3 Computation of the modulation equations. 

We now compute the modulation parameters using the sets of orthogonality con- 
ditions and the e equation (|3.1ip . (|3.12p . This is done in detail in Appendix A, and 
the bounds (IH?^ . (CmD follow. 
This concludes the proof of Lemma 13.11 

3.2. Refined energy identity. Our aim in this subsection is to derive a general 
refined mixed energy/Morawetz type estimate which will allow us to derive a Lya- 
pounov function for critical mass blow up solutions. We shall work within the 
following general framework. We let u be a solution to (jl.ip on [tQ,0) and w be an 
approximate solution to (jl.ip : 

idtw + Aw + k{x)\w\'^w = il^, (3.31) 

with the a priori bounds 

< 1, \\^w\\l2<^, II-wII 3 < (3.32) 

We then decompose u = w + u so that u satisfies: 

idtu + An + k{x){\u\'^u — \w\'^w) = —ip (3.33) 
and assume the a priori bounds on il^: 

l|Vn||i2<A, ||n||L2<A2 (3.34) 
and on the geometrical parameters: 

|AAt + 6|<A^ 6~A, |Aat|<A, |6t| < 1 (3.35) 
for some nonnegative parameters < A, 6 ^ 1 and a G M^. 

We let A > be a large enough constant which will be chosen later and let (j) : 
— )• M be a smooth radially symmetric cut off function with 



Let 



F(n) = ^|n|^, /(n) = n|n|2 so that F\u) ■ h = Re{f{u)h). 



We claim the following generalized energy estimate on the linearized flow (|3.33p : 



14n 



The standard orthogonality condition in p.30|l on £2 is (e2,Q) = because Q is the bound 
state L-Q = 0. In |26) . it is shown how to replace this by (£2,A^(5) = 0, see proof of Lemma 
3, and the same proof applies verbatim with the choice (e2,p) = using the key non degeneracy 
dill). 

^^corresponding to the asymptotic regime near the singularity, see (|4.60|l 
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Lemma 3.3 (Generalized energy estimate). Let 



+ 



1 
2 

1 b 

2 A 



ivnr + 



1 



2 



u\ 

X — a 
AX 



k{x) [F{w + u)- F{w) - F'{w) ■ u] 



Vuu 



(3.37) 



then there holds: 



dl 



— = {^J k{x)w'^u'^ ] — 5ft ( / k{x)wt{2\u\'^w + u'^w) 



4^2 



A 



(3.38) 



2 . I X — a \ |tt'2 



/ A2 



b ^ f I , I X — a \ , , w„,-|2 ~2— \ 



+ An^'^^ ^a 



fc(a;)(2|{t| w + u w) ■ Vw 



All! - ^ + k{x){2\w\'^^ - v?iIj) + i\AV(j) ( ^^^-p^ 
A Ay AX 



. h ^ f X — a . , 



n)+0(A2||V^||2, 



1 



I 1 1 2 ~ 



Remark 3.4. The virtue of (|3.38p is to keep track of the quadratic terms in u. In the 
various situations we will encounter, the key will be to prove that the corresponding 
quadratic forms both in the boundary term in time and in the RHS of ()3.38p are 
small deformations of the linearized energy and hence are definite positive. This 
will generate a control of the form: 



L2 



r,ll2 



\U 



L2 



A2 



b 

^1? 



|x-a|<A 



|Vn|2 + 



lL2 



A2 



+ l.o.t. 



(3.39) 



The two keys are first the sign of the above RHS and second the fact that this 
is a useful information once the blow up speed is known, that is 6 ~ A. The 
Lyapounov property (j3.38p will be the key to estimate the solution backwards from 
the singularity, which will be used both for the proof of existence and uniqueness. 



Proof of Lemma 13.31 

step 1 Algebraic derivation of the energetic part. 
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We compute from (|3.33p : 



1 f \u\ 



2 / A2 



k{x) [{F{u) - F{w) - F'{w) ■ u)] \ (3.40) 



-3f? ( dtu, Au-^u + k{x){f{u) - f{w))\ - ^ / |u|2 



3f? dtw, k{x){f{u + w)- f{w) - f'{w) ■ u) 



(^V, ^u-^u + k{x){f{u) - f{w))^ - ^9 (fc(x)(/(n) - /H), n) 



- ^ y \u\^ - (^dtw, k{x){f{u + w)- f{w) - f'{w) ■ u)) 

= 'is Au — + k{x){2\w\'^u + uw"^)^ — J k{x)u'w^ 

~ ^ j - ^ {dtw, k{x){wifi + 2u;|n|2)) 

+ 9 ~ 3^^' ^(^)(/(^ + u) — f{w) — f'{w) ■ tt)^ — (j)tw, /c(x)u|ttp^ 



where we used that f'{w) ■ u = 2\w\'^u + w^u. We first estimate from (I3.35P : 



- ^ / I^P = ^ / - ^(AA, + b)\u\l, = ^ l\u{' + O(||n||^0. (3.41) 



It remains to estimate the last line in the RHS (|3.4Up . For the quadratic and higher 
terms, we estimate using the a priori bounds (|3.32p . (|3.34p : 



1 



V' - ^u, k{x){f{w + u) - f{w) - f'{w) ■ u) 



^ [il) — ^^^) k[x){v?w + 2|M|2r(; + lupn) 



< 



< 



L2(M2)||n||i6(||w||L6 + II-uIIlO + (1 + IklliOII'" 
1 ii2/3ii ii4/3 ii2 



,~,I|3 
L6 



< A^IIV-lli^ + Nl^i. 



(3.42) 



For the cubic terms hitting wt^ we replace wt using (|3.3ip . integrate by parts and 
then rely on ()3.32p to estimate: 



/ 



k{x)wt\u\'^u 



^ \\w\\h3/2\\\u\^u\\h^/2^^^2) + \\w\\\(i\\u\\\6 + ||V'IIl2(k2)||'u|||6 



< ^ ll.-,llV2||~||5/2 , 1 II-II2 



A3/2 



\u\ 



L2 



+ + IIV'llL2||n||^l||u||i2 



L-;l|2 



(3.43) 
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Injecting (|3.4ip . (|3.42p . (|3.43p into (j3.40p yields the preliminary computation: 



d r 1 r^,,, 1 f 



k{x){F{u) - F{w) - F'{w) ■ u) 



+ 9 



lb — 
^V' - + k(x)(2\w\^i; - w'^tp) 



u 



+ o{x^\mi, + -\\u\\i, + \\u^^' 



A2' 



lL2 



(3.44) 



step 2 Algebraic derivation of the localized virial part. 
Let 



^7/ N b , f X — a 



then 



11 (b^ 
2dt\X' 

1, 



(^) ^^^) 



+ • Vu 



(3.45) 



Using ()3.35p . we estimate in brute force: 

1 



< -^{\X\ + b^\ + \XXt + b\) + ^{\Xat\ +b)<l 



from which: 



< j\\u\\l2\\Vu\\l^ = O (^\\u\\l2 + Mil ] . (3.46) 



The second term in ()3.45p corresponds to the localized Morawetz multiplier, and we 
get from (|3.33p and integration by parts: 



3? I / idtu 
b 



-A(/>n + V(j)-Vu 



(3.47) 



( / V^'^ ( ^ ^ (Vti, W) ] - f / A^' 



' 7 AA^X^ 
k{x){\u\\ - \w\^w) ■ 



X — a 
AX 



I b ^ f f . , f X — a 



k{x){\u\'^u — \w\'^w)u 
1 b 



X — a \ — 
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We now expand the nonlinear terms and estimate the cubic and higher terms: 



A 



X — a 



< 



2A2^ 

ni„--,ii3 



X — a 



k{x){2\u\ w + u w + \u\ u) ■ Vu 



k{x){2\u\'''w + u^w + \u\'''u)u 



L~.|2- 



1 



^IIl6(k2) + ||^t|lL6(]R2)||?w||L6(lR2))||Vu||i2(R2) + y ( 1 1 U 1 1 ^4 (IR2) + ||n 



,~l|3 



< ll~l|4 II ||l/3 II ||2/3 

\m\m(M.2) + |F||^l(lR2)||'U|l2,2(K2)||W^|li2(]R2)||W^||^l(K2) 



< lU 



~||2 I ||~||3/2 



1 3/2 



1 1/2 



1 1/2 



where we have used (|3.34p and ()3.32p . 

The remaining quadratic terms in (|3.47p are integrated by parts: 



A 



X — a 



AX 



.h , f X — a\ „, . h . , ( X — a . , 



(3.48) 



A 



( 



^ " ^ A;(x)(2|u;pn + w^u) • Vn ) 
/i;(x)(2|ty|^-Lt + w^u)u 



\ A\ 

X — a 



A 



2A 



^A 

X — a 
AX 
X — a 
AX 



k{x){2\u\ w + u w) ■ Vw 



The last term gains an extra A smallness: 

'avJ'-"" 



X 



I^IIl4||^Iil4 



V AX 

1 



Vk{x) 2|np|tt;p + M^tM^ 



Vk{x) \2\u\'^\w\'^ + ifw'^ 



< WuW'i^W in 111 A < — II Vn||^2 1111112,2 = O I \\u\ 



2 -L 



U 



L2 



(3.49) 



(3.50) 



Injecting (|3.48p . (|3.49ll . p.50p into ()3.47p yields after a further integration by parts: 



+ 



5? 

b_ 
A2 

b 



idtu 



^A0u + V0 • Vn 



A 



X — a 
AX 



{Vu,Vu) 



1 b 



4A2A4 



A' 



X — a \ , ,9 



AX 



AV(j) 



/ X — a ^ ^ ^ \,'r)l ~|2 



+ 9 



i-AV(j) 
A 



V AA 

X — a 



AX 



A;(x)(2|tt| w + u w) ■ Vw 



^, ■ b ^ , f X — a , , 



u] +0\\\u 



~l|2 



+ 



A2 



Together with (|3.44p - (j3.46p . this yields p.38p . This concludes the proof of Lemma 
13:31 
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3.3. Backwards propagation of smallness. The first application of the energy 
estimate (|3.38p is a bootstrap control on critical mass solutions to (jl.ip . More 
precisely, let uhe a solution to (jl.ip defined on [tQ,0). Let to < ti < and assume 
that u admits on [to^^i] a geometrical decomposition on the form: 

where e satisfies the orthogonality conditions (|3.4p - (j3.8p and +|7^(t)| ^ 1. 

Let 

u{t,x) = -^-^^,e ( t, ^^^) e^-W. (3.51) 
Assume that the energy £^0 satisfies the lower bound and define Cq as: 



WyQf 



L2 



(3.52) 



^ 8Eo + fVk^O){y,y)Q^- 
We claim the following backwards propagation estimates: 

Lemma 3.5 (Backwards propagation of smallness). Assuming that there holds for 
some ii < close enough to 0: 

\\\uU2-\\Q\\L2\<X\t,), (3.53) 



|Vn(ti)||i2 + 



+ 



ih) <A(ti), 



A(ti) + 



ti_ 



<A^(ii), 



A(ti 



1 



(3.54) 
<X\t,). (3.55) 



Then there exists a backwards time to depending only on Cq such that\/t € [to,ii], 

m 



iVtt 



L2 + 



A2(t) 



?^<l|Vn(tOlli. + ^^ + A«(t) 



a(^)-^ 



{t)\<\\ti 



L2 

<A^W, 
< A^W, 



A2(f 



(0 



\{t). 



(3.56) 
(3.57) 
(3.58) 
(3.59) 



Proof of Lemma 13.51 



Let us consider from the continuity u G C([to, ^i], ^^) a backwards time to such 
that Vt G [to, ti], 

||ti(t)||^2 < K\\t), \\u{t)\\m < K\{t), (3.60) 



m 



+ 



a(t) 



A(t) 



< KX{t), 



< i^A3(t), 



Kt) 



1 



A(t) Co 



<K\^{t). (3.61) 



for some large enough universal constant K > Q. Then we claim that p.56p . 
(j3.57p . (j3.59p hold on [to,ti] hence improving the bounds (j3.60p . (j3.6ip on [to>ii] 
for to = ^o(Co) small enough independent of ti. 



step 1 Monotonicity of the norm. 



-'^''which will be proved in (|4.17|l 
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Let us apply Lemma [3.31 with: 
w{t, x) = Q{t, x) 



1 



1 



■Qr{t) 



x-a{t) \ .^(4) 

m ' 



(3.62) 



Let X be given by (|3.37p . we claim that ()3.38p implies the following coercivity 
property: 

f > A / \n\^ + 0{K'X' + \\u\\U (3.63) 



Assume (|3.63p . By Sobolev embedding and the smallness of e, we have the rough 
upper bound on Z: 

(3.64) 



Now after renormalization, the proximity of Qp to Q ensures that the quadratic 
part of X is a small deformation of the linearized energy and hence the coercitivity 
property ()3.30p with the choice of orthogonality conditions on e ensures : 



1 
2 

2 A' 



2 



A2 
X — a 
AX 



Vuu 



(3.65) 



= ^ [{L+ei,ei) + {1.82,82) + o{\\e\\l,)]>^ J \\e\\l, - {e,,Qf 

We next estimate from the conservation of the norm ()3.22p . ()3.53p . (|3.6ip and 
(13:23]) : 



<\\8\\l,+K'X\t) 



Qv)\ < Mlh + X\t) + + J - J Q' 
from which 

iei,Q)'<oi\\8\\l,)+K'X'it). (3.66) 
We then integrate (|3.63p in time and inject (j3.64p . p.65p . (j3.66p to conclude: 



iVn 



L2 



+ 



^ < l|Vn(ti)||i.+ 



Hti)\\ 



L2 



+ K^X^{t) + 



A2(i) 
u{ti)\\] 



< l|Vn(ti)||i. + ^^ + /^^A6(t) + ^ 



|n(r)||^,+i^^A^(r))cir 



|Vn(r)||i.+ 



A2(r) 



for to = to (Co) small enough, and (I3.56P follows from Gronwall's lemma. It implies 
in particular together with (|3.54p : 



\Vu{t)\\l,+ 



(3.67) 



and closes the bootstrap of (|3.60p . 

step 2 Integration of the law for the parameters. 

We now integrate the law for the parameters. Indeed, (|3.20p . ()3.6ip and (|3.67p 
imply like for the proof of ()4.80p : 



\bs + b - do{a, a)\ + 



A. 



+ b 



+ |/3, + 6/3-co(a)A-/33A3| + 



T 



2(3 



< X^ + K^X"- 

< X\ (3.68) 
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One should be careful when reintegrating the above system which is stable thanks 
to (j3.68p and the sign do (a, a) < provided by (I3.18P : 

b\ bs + b"^ — do{a,a) b ( Xg 
A, 



A 



4(a,a) < ^3 



A 



and hence: 



1 



^-^(^)<^-^(si) + A2(.)<a2(.) 



(3.69) 



where we used (|3.55|) . We now write down the conservation of energy at t and add 
the conservation of L? norm. We get using verbatim the same algebra like the one 
which led to (f3:28|) and using (jS^GT]) : 



1 



V^miy,y)Q^] +kia 



and thus from the choice of Co (|3.52p and (|3.53p : 



+ 



a 



1 



1 



A2 ' A2 ~ C72 A2 ^ ~ Co A ~ 



where we used (|3.69p in the last step. This together with (j3.69p again yields: 



<A^ 



b 1 
A " Cb 

We eventually have from (|3.68p for the scaling parameter: 

b 



(3.70) 



and hence using (I3.55p . (|3.70p : 



m + TT 



< 



+ 



b 1 



dT + xHt)<X^it). 



This concludes the proof of (I3.57p . (|3.58p . (|3.59p assuming (|3.63p . 
step 3 Coercivity of the quadratic form in the RHS of (|3.38p . 



We now turn to the proof of the Lyapounov property (|3.63p and start with com- 
puting more explicitly the quadratic terms in the RHS of (|3.38p for w = Q given by 



/C(n) 



we claim: 



b 

A2 



A2 



X — a 
AX 



k(x)wt(2\u\'^w + u'^w) 
1 



4^2 



X — a 



k{x){2\u\''^w + u^w) ■ Vw 



/C(n) > 



A3 



Ve|"e~^ + / |e|" ) +0{K^X^) 



A' 



X — a 



AX / A2 



(3.71) 



(3.72) 
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for some universal constant c > 0. Indeed, first observe from (j3.18p . p.20p . (|3.60p 
and (13:61]) : 



\Mod{t)\ + \do{a,a)\ + \di{a,a)\ < K^X\ 
We then first compute from ()3.62p : 



(3.73) 



Qt 



(Xat-2l3)Vk(a) 2/3VA;(a) XXt + b b i X^^t 



2k{a)X 



2k{a)X 



A2 ' A2 ' A2 



A2 



+ 



+ 



AAf + b / X — a\ b fx — a\ Xat — 2/3 2/3 \ 



A2 
1 



A 



A2 \ A 



A2 



A[A;(a)]5 



X — a{t) ' 



fc(a(t))V2A(t) * 5P V Kt) 



i b \ ~ b / X — a 

A^+A^r^+Ahr 



VQ + ( ye 



where we used ()3.73p . ()3.6ip in the last step. We may thus use the exponential 
decay of Q-p to conclude after renormalization: 



K( I k{x)Qt{2\u\^Q + u^Q) ] = { I k{x)Q{2\u\^Q + u^Q) 



A2^ 
1 



;?fi[ I k{x){2\ufQ + u^Q)Q 
+ ^^0{KX\\e\\l,). 



b 



X — a 



k{x){2\u\'^Q + u^Q)-VQ 



We now inject this estimate into ()3.7ip and rewrite the result in renormalized vari- 
ables: 



/C(n) 



- Jk{a + Xy){{\Qv\^ + 2^^)sl + 4SGei£2 + {\Qv\^ + 2Q^)el) ^ 



2 J. (]r 



+ — 



AX 



y] k{a + Xy){2\e\^Q'p + e^Qv)-^Qv 



+ -0{KX\\e\\l,). 



From the proximity of Q-p to Q, the above quadratic form is for A large enough a 
small deformation of the localized in A linearized energy, and hence (|3.30p . (|3.6ip 
and our choice of orthogonality conditions ensure for A large enough: 



/C(n) > ^ 



|Ve|2e VA + / |e|2 _ (^^^q)^ 



- A3 



iVere ^ + \e 



+0{K^X^^ 



where we used (|3.66p in the last step, this is (j3.72p . 

step 4 Control of the remainder terms in the RHS of (j3.38p . 
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It remains to control the ■0 terms in (j3.38p . According to (I3.3ip . p.62p and the 
construction of Q-p, we have: 



1 



A3A;(a)i/2 



a,9 



i{bs + b')dbQv + M y + M ^^aQp 



(3.74) 



A. 



A 



i [-j- - 2/3j Aa„Qp + i(/3, + 6/3 - Aco(a) - /SgA^ - PAX')dpQv -i[^ + h]AQv 
^ - 2/3) . {vQv + \^^Qv ]+{%- \P?)Qv + i^v 



where ipp is the remainder in the construction of Q-p: 

Let us start with a rough bound on ip: from (I3.20p . (j3.60p . (|3.6ip . there holds for 
i = 0,1,2: 



e [|Mod(t)| + |a|2 + A^] 



1 



< 



\x — a(t)\ 



and thus: 

This yields in particular the bounds: 

>^'mh<\\e\\h+K^X', 



(3.75) 



(3.76) 



(3.77) 



< 



< 



9 j (^ijAV^ 

||VV'||L2||^2||i2 + 



X — a\ . b ^ , f X — a , , . 



L2 



if 



L2<^[A||e||L2+A^] \\eh2 



IL2 



A3 



+ if4A5. 



(3.78) 



The rough bound (|3.75p is not enough to take care of the remainder term in (|3.38p 
for which we need a further cancellation. From (|3.74p and the construction of Q-p, 
we have ip = ipi + ■ip2 with 



1 



—I 



A. 



1^' 
4 



{bs + b^Y-^Q - {(3s + 6/3 - Aco(a) - /SgA^ - (3^X%Q 



a. 



^ +6 1 AQ-i^:j_2/3) VQ + (7.-|/3|')g 



m 



(3.79) 



and for i = 0, 1, 2: 

1 \x-a(,t)\ 

|VV2(x)| < T^e— 



X3+i 



^[\r\Mod{t)+KH^]< ^^^^ 



- MO [A^l 



e||i2 +if2^5i 
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The contribution of this remainder term 'tp2 is estimated in brute force: 

'^2 , ,2 , 2- 



AV'2 + k{x){2\w\^'il:2 - w V'2) ) u 



< 



I|AV'2||l2 + 



A2 



2 



We now observe that the ipi term hes after renormalization in the generahzed nuU 
space of L and thus a 0{V) factor is gained using the choice of orthogonahty con- 
ditions on e. In other words, after renormahzation: 

A^i — ^ + k{x){2\w\'^'ipi — w'^ipij u 



< 



< 



< 



\Mod{t)\ + \a\ 
A4 



(e2,i-(|2/PQ))| + \{e2,L.{yQ))\ + L+(VQ))| 



+ |(£2,i-Q)|+0(|^|||e||LO] +T4 



A. 



2\4 



\\\e\\L2+ K^\ 



+ h 



2 \5 



|(ei,L+(AQ))| 



\V\\\e\\j2 + 



\\\e\\L2+ K^\ 



2\4\ 



e\\L2 +K^X 



A3 



(3.80) 



where we have used in particular the estimate (|3.2ip and the conservation of mass 
(13:66]) to bound the term + b\\{ei, L+ ( AQ) ) | . 

We now inject (IHTTD . pTTHjl . (prHUjl and (IHT^ into (IH:H5D and (pmjl fohows. 
This concludes the proof of Lemma 13.51 



3.4. Existence of critical mass blov^ up solutions. A first consequence of the 
backwards propagation estimates of Lemma 13.51 is that one can integrate the flow 
backwards from the singularity to prove the existence of critical mass finite time 
blow up solutions. This approach is similar to the one in |24] . [19| . 1221 . 1151 



Proposition 3.6 (Existence of critical mass blow up solutions). Let 

70 GM, Eo>^J V^m{y,y)Q\ 

and Co given by (|1.13p . then there exists tQ < and a solution G C([to,0),// 
to (jl.ip which blows up at T = with 

E{uc) = Eq and \\uc\\l2 = WQWl'^- 
Moreover, the solution admits on [to,0) a geometrical decomposition: 



Qc + Uc (3.81) 

where Ec satisfies the orthogonality conditions ()3.4p - ()3.8p . and there holds the bounds: 

\ur\\m<\l ll^iclL,3 < a|, (3.82) 

(3.83) 



\Uc\\l^ < Ac, 



t 



br. 1 



K + -^ = 0{\i), ^- =0{Xi), 

Go Ac Oo 



l«c| + |/3c| < Xl 7c 



fi2 



+ 70 + 0(Ac 



(3.84) 
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Proof of Proposition 13.61 



step 1 Backwards uniform bounds. 



Let a sequence t„ — )■ and u„ be the solution to (|l.ip with initial data at t = tr, 
given by: 



with Vn = {bn{tn)An{tn),Pn{tn),an{tn)) and: 



^0 

We have by ([225]): 



, an{tn) = Pn{tn) = 0, 7n(in) = 70 



r2 



(3.85) 



(3.86) 



|n„(t„)||i. = / Q' + 0{tl). 



and Un{tn) = by construction. Hence Un satisfies at ti = tn the assumptions of 
Lemma [3.5[ and thus we can find a time to independent of n such that Vt € [to,tn), 
Un admits a geometrical decomposition 

1 1 



Un{t,x) 



[A:(a„(t))]U.(t)^^"(n'' Kit) 



with uniform bounds in n: 



iVu 



-it)-- 

An Co 



t 



A^(t) 



< A^W 



(3.87) 



An(0 + 



<XUt), |«n(t)| + |/3n(t)|<A2(t). (3.88) 



From standard Strichartz bounds, this implies the uniform bound: 

3 

\\Un\\ ,, , 3 <, A„(t) 2 , 



(3.89) 



which we prove in step 2. 

The compactness of Un{to) is now a consequence of a standard localization 
procedure. Indeed, let a cut off function x(^) = for |x| < 1 and x{^) = 1 ^or 
\x\ > 2 and xr{x) = x(^), then 



/m / VxR • Vn 



<1 



where we used (|3.87p . (|3.88p and (|3.89p . Integrating this backwards from ti to to 
using ()3.85p . (|3.86p and ()3.89p yields up to a subsequence: 

Unito) —J" Uc(to) in as n — )• +oo. 

Let then Uc be the solution to (II. ip with initial data tic(to), then the continuity 
of the flow ensures: Vt G [to,0), 

Un{t) — ?• Uc{t) in 
and tic admits a geometrical decomposition (|3.8ip with 
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see [26], [20], [15] for related statements. The bound (|3.82p and the estimate on 
the geometrical parameters (|3.83p follow by passing to the limit in (|3.87p . (j3.88p . 
This implies in particular that Uc blows up at t = 0. Eventually, the conservation 
of the norm ensures: 



|iic||i2 = lim ||u„(t„)||i2 = IIQII 



L2- 



Similarly, 

E{uc{t)) 



2X1 



k{x)\uc{t,x)\' 

2 



\yQ\h 

8 A2 



-^J kiX,y + a,)\Qpf + OiX, 
l-^Jv^k{0){y,y)Q^ + o{l)^Eo as t ^ 



and hence the conservation of energy and the choice of Cq (|1.13p ensure: 



E{uc) = Eo. 



Eventually, we derive from (|3.87p . (|3.88p the rough bound 

\ifn)s\<Xl 

which implies using (|3.88p : 

1 



d ( cl 



1 
AH 



(in 



f<2\2 



t2 



XI 



2\2 



Cn A: 



t2 



< 1. 



We then integrate in time using (|3.86p to conclude: 

(J2 

7n(t) + ^ = 7o + 0(i) 



and hence the limit in the phase parameter (|3.84p follows by taking the limit 
n — )• +00. 



step 2 H2 bound. 



It remains to prove the ffz bound (j3.89p which implies in particular the bound 

3 

(j3.82p by passing to the weak limit. 
Un satisfies: 



idtUn + An„ = -iln - k{x)Un\Ur, 



with 



idtQn + ^Qn+k{x)Qn\Qn? = V'n, En = k{x){Qn+Un)\Qn+Un\^ -k{x)Qn\Qn? -k{x) 

Hence from standard Strichartz bounds and the smoothing effect of the linear 
Schrodinger flow, there holds using u{tn) = 0: 

llV^UnlU- ,L2 ^ llV^V'nlU/a 4/3 + 11(1 + |x|2)F„||^2 „i 



+ ||V2(n„|tt„|2) 



(3.90) 



The error if^n is estimated from ()3.74p which yields a bound: 



VtV'nll^l < ^{\Modn{t)\ + lanP + A^) < A„ 
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with Modn{t) corresponding to the modulation equation of Vn{t) hke for p.20p . 
and where we used from IK20\i . dSSZD, ([3:881) : 

Modnit) < Xt 

Hence: 

||viV'n||^4/3 ^4/3 <Ai (3.91) 

The Fn term is a local term in y which contains linear and quadratic terms in u, 
hence from (|3.87p : 

1111 ^3 

||(l + |x|^)F„||j:^i < -rT\\Un\\L^ + l-2\\Un\\m+l-2\\Un\\% + l-\\u^^^ < A„+A^, || V 2 M„ || ^2 

A„ A„ A„ An 

and thus: 

||(l + |x|2)Fn|lL2 < Xi + XlWVhnh?^^ r2. (3.92) 

Eventually, the nonlinear term is estimated from Sobolev embeddings and standard 
nonlinear estimates in Besov spaces: 

||vi(^2„|n„|2)||^4/3 < \\V'/^Unh4unfHi < ||V3/2^i„||^2At 
Injecting this together with (f3M]l . (fSlMjl into (fSlQO]) yields: 

and ([3:89]) follows. 

This concludes the proof of Proposition 13.61 



4. Critical mass blow up solutions have conformal speed 

We now turn to the proof of uniqueness of the critical mass blow up solutions at 
a nondegenerate point. We let k satisfy Assumption (H) and let u(t) G be a 
solution to (|l.ip with critical mass 

\Ml^ = IIQIIl2 

and which blows up at T = 0. We let Eq = E{uo). 

We start with slightly revisiting the variational arguments of [25\ to show that the 
solution admits near blow up time a decomposition: 

uit x) = - -Q (t ^^^^ e^7W + ^ 

with the following a priori bounds: 

||n(t)||^i <C7(no), (4.1) 
a{t) ^ a* as t with k{a*) = l, (4.2) 

A(t) < C{uo)\t\ ie \Vu{t)\L2 > ^1^. (4.3) 

In particular, the solution will be in the regime described in section 13.11 
The main difficulty at this stage is that there holds no a priori upper bound on the 
blow up rate, and the remaining subsections are devoted to proving that u must in 
fact blow up with the conformal speed: 

l|Vn(t)|U2~^ (4.4) 
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which is a key step towards classification. The key will be to exhibit a new Lya- 
pomiov type rigidity property adapted to the critical mass setting -Proposition 14.31 - 
which will allow us to improve the energy bound (|4.ip for a dispersive bound: 

\\uit)\\m ^0 as t 0. 

4.1. Variational estimates and convergence of the concentration point. 

We start with revisiting the variational estimates in [25j in order to derive the 
convergence of the concentration point and the lower bound on the blow up speed. 

Lemma 4.1 (Variational control of minimal mass blow up solutions). Let u{t) be a 
critical mass solution to (jl.ip which blows up atT = 0. Then for t < close enough 
to 0, u{t) admits a geometrical decomposition 

u{<..) = ^«3 + .)(<,^).-'<". (4^5) 

for some parameters {X{t), a{t),^{t)) G M!i_ x x M with: 

1. Uniform bound on the decomposition: 

|l-fc(a(t))|V2 + ||e(t)||^i(jj2^ <A(t)^0 as t^O. (4.6) 

2. Convergence of the concentration point: 

a{t) a* with k{a*) = 1. (4.7) 

3. Lower bound on the blow up rate: 

X{t)<C{uo)\t\. (4.8) 

Proof of Lemma \4-l\ 



The proof is standard and relies on the orbital stability of the ground state soli- 
tary wave. We will briefly sketch the argument and the computations which are a 
simplified version of the proof of Lemma 13.11 



step 1 Nonlinear decomposition of the flow. 



Let 

vo{t,x) = Xo{t)u{t,Xo{t)x) with Xo{t) - '^^'^'^ 



|Vn(t)|i2' 
then 

\voit)\L2 = \u{t)\L2 = \Q\l2, \Vvo{t)\L2 = \VQ\l2 
and the conservation of energy and Assumption (H) imply 

I J |Vi;o|'-^ j \vo\^ <\j \'^vo?-\ j k{X^{t)y)\vo\^ = Xlit)Eo ^ as t ^ 0. 



Hence from a standard concentration compactness argument, see |18| and Lemma 
1 in [26], there exist (xo(t), 7o(t)) G x M such that 

vo{t, ■ + xo(t))e*^oW ^ Q in as t ^ 0. 

Hence u{t) admits near blow up time a decomposition 

x-a{t) \ ^i^(j) 



«(M) = ^(0 + .)(^«,^j.-™ (4,9) 
with 

\l-k{a{t))\ + \\e{t)\\Hi ^0 as t ^ 0. (4.10) 
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Using the implicit function theorem, the uniqueness of the decomposition (|4.10p can 
be ensured through a suitable choice of orthogonality condition4_|. We then set the 
orthogonality conditions on e to be: 

(ei,|2/|2Q) = 0, (£i,yQ) = 0, (£2,p) = (4.11) 

where we wrote 

e = ei + ie2. (4.12) 
Let V = Q + then v satisfies an equation similar to ()2.2p : 



idsV + — V + k{\{t)y + a{t))v\v\'^ = i-^Av + i— ^ • Vv + 

A A 

with 7s = 7s — 1. We next write down the equation satisfied by e and compute 
the geometrical parameters using the orthogonality conditions ()4.1ip . see [26j for 
very closely related computations or Appendix A for more details in the setting of 
Lemma 13.11 A simple computation left to the reader leads to the regularity in 
time of the geometrical parameters together with the bound: 



+ 



A. 
A 

step 2 Expansion of the conservation laws 



A 



<\e\Hi + \l-k{am+Xit). (4.13) 



We now expand the conservation laws in the e variables. For the norm, we 
have from the critical mass assumption and (|4.9p : 

j \u{t,x)\^dx = j \v{s,y)\'^dy = j 

and thus: 

2(ei,Q) + j kl' = 0. (4.14) 
For the energy, we have by rescaling: 

A^i^o = \ J \yv{s, y)\^dy - i / HXy + a)\vis, y)\Uy = E{v). (4.15) 
We compute on Q which has zero energy: 

E{Q) = \ j fc(Ay + a)g' 



where used in the last step the boundedness of k' s derivatives and the radiallity of 
Q. 

We then expand ()4.15p and inject the conservation of norm of ()4.14p to get 
after a classical computation: 

2\^Eq = ^j{l-k{a))Q^ + {L+ei,ei) + (1^82,82) 

+ OiX^ + {X + \l-k{a)\)\\8\\Hi + \\£fHi). (4.16) 



"'^'''Note that the profile Q corresponds to a very rough approximate solution to 1)2. 4|l . but at 
this stage, we do not know yet that a is stabilizing around a* and hence we cannot introduce the 
refined profiles Q-p, V — {b,\,l3,a) which make sense for a close to a* only. 
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The coercivity property of Lemma 13.21 together with our choice of orthogonahty 
conditions (|4.1ip and the degeneracy inherited from (|4.14p : 

\{euQ)\' = o{\\e\\l,) 

yields: 

(L+ei,ei) + (L_e2,e2) > y||e||^i- 
Injecting this into ()4.16p together with the a priori smahness ()4.10p yields: 

(l-fe(a)) + ||e||^i <A2, 

and (14.61) follows. 



step 3 Convergence on the concentration point and upper bound on the blow up 
rate. 

We now prove (|4.7p . ()4.8p which are a straightforward consequence of ()4.6p and 
()4.13p . Indeed, they first imply 



\at\ = J 



A 



/ ||e||^^i +X + \l-k{a)\ ^ ^ , , /"° , , , 
< ^^-^ ^— < 1 and thus / \at\dt < +oo. 



'-1 

Hence a{t) — )■ q* as t — t- and k{a*) = 1 from (|4.6p . This proves (|4.7p . Next, from 
gSD and (|¥T3D . 

A, 



n4 



A 



, Ikllffi + A + |1 — k(a)\ . , , , X / ^ / I I 
^\U1H _\ L^<1 and thus X{t)<\t\ 



by integration from t to and the blow up assumption at t = 0, i.e. A(0) = 0. 
This concludes the proof of Lemma 14.11 

4.2. Strict lower bound on the energy. From Lemma 14. 1| we know that the 
center of mass of u must stabilize around a point a* where k reaches its maximum. 
Without loss of generality with respect to the assumptions of Theorem II. 3| we may 
assume that 

a* = 0, A:(0) = 1, V/c(0) = 0, V^k{0) < 0. 

Hence the solution is from (|4.6p on [tg, 0) for to close enough to blow up time in the 
regime described in section [3Tl and we may now sharpen the decomposition (|4.9p 
by introducing the modified Qp profiles given by ()2.2ip which by construction are 
a small deformation of Q. We thus decompose the solution for t close to enough to 
as 

where e satisfies the orthogonality conditions ()3.4p - ()3.8p . Moreover, from (|4.8p . the 
rescaled time ()3.10p is a global time: 

= / TTTT as t -)• 0, 

J to ^ iV 

and thus the solution satisfies the bounds (|3.17p . (j3.20p of Lemma 13.11 Vs £ 

[so,+oo), 

52 + |^|2 + |^|2 ^ ||e||2^^^^,^ <x^(^Eo + ^J V^k{0){y, y)Q^^ + 0{V^), 
Modit) <V^ + V^\\£\\l2(w?) + \\e\\l, + \\e\\% + {h-^^V'' + V{a^ + p^). 
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We now claim that this implies in particular the strict lower bound on the energy: 
Lemma 4.2 (Strict lower bound on the energy). There holds: 

Eo>-\j V^k{0)iy,y)Q\ (4.17) 

Proof of Lemma 14.21 

By contradiction, if Eq + If V'^k{0){y,y)Q'^ < 0, then (l3T7l) yields: 

62 + |«|2 + | ,2^n ||2 <x4 



which together with (I3.20p . implies: 



Dividing by A^, we get: 



^ = 0(1) ie ^ = on) and hence ln(A(t)) = 0(1) as t ^ 0. 

But this contradicts the fact that A(t) — )• as t — )• 0, and ()4.17p is proved. This 
concludes the proof of Lemma 14.21 



4.3. The localized virial identity. We now turn to the heart of the proof which 
is an improved local bound for e locally on the singularity. Our aim is to strictly 
gain on the global in space energy bound ()3.17p . Let us stress that the fact that 
there holds no a priori upper bound on the blow up rate means that one cannot rule 
out a priori a regime for which 

6 < A, 

where we recall that 5 ~ A is the expected pseudo-conformal regime. This makes 
the refined energy identity (j3.38p useless at this stage because of the presence of 
the b factor in front of the quadratic term of the RHS of (j3.39p which may lead this 
gain to degenerate. 

We claim that a Morawetz type computation in the spirit of the local virial estimate 
initiated by Martel and Merle for the study of the mass critical KdV [20j allows 
one to obtain a space time bound for weighted norms centered on the singularity 
which improves the energy bound (j3.17p . Let us however insist onto the fact that 
a very specific algebra is at hand here which relies on the critical mass assumption 
and avoids both the standard troubles when running two dimensional Morawetz 
estimates for Schrodinger, and the delicate study of explicit quadratic forms like in 
[20| . [26j, [9], the only required coercivity property here being the energetic lower 

bound daini). 

Proposition 4.3 (Local virial control). Let (j) be given by (|3.36p . There exist uni- 
versal constants c, c > and a large enough constant ^ > such that for t close 
enough to 0, we have: 
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The strength of the estimate (j4.18p is that the smaller the A -i.e. the faster 
the blow up speed-, the better is the estimate. Indeed, the terms involved in the 
boundary term in time are uniformly bounded from (|3.17p : 



\yQ\ 



L2 



1 



2A 



+ —9 / [^]Vee 



y 



A 



<\b\^mL<i 
~ A A ~ ' 



(4.19) 



while from ()3.17p again and the finite time blow assumption: 

{Mm+r' + m')ds< / x'{s)ds< / dt<i. 

Hence, integrating (|4.18p between sq and +oo, and using (|4.19p (|4.20p yields: 



(4.20) 



so 



-OO -j^ 

A 



On the other hand, 



+ 00 



X{s)ds 



to 



dr 

W) 



+00 



(4.21) 



(4.22) 



from ()4.8p and hence ()4.2ip is indeed a strict gain with respect to the energy bound 
Proof of Proposition 14.31 

The proof relies on an algebraic computation and in particular the specific struc- 
ture of the quadratic terms in £ appearing in the RHS of ()4.18p . 

step 1 Computation of {j) ^■ 

We have: 

^\ _ bs + b^ b fx, 

aA a aVa 

which together with the law ([SlMP of 60, and (jSTT]) . ([3:20]) yields: 



\y\ 



■Q,AQ 



V2fc(0)(a,a) 



+ 



(i?i(e),AS)-(i?2(e),Ae) 



+ O {V^ + V\\£\\l2 + ||e||^i + lap + . 
Using the definition ()3.15p of Ri and ()3.16p of i?2, we obtain: 



^)^(M!,,.,)4[_Z!M(^ 

- (3Sef + 2eeie2 + Se^, AS) - (396^ + 2Seie2 + Oe?, AG)] ) 
+ 0{V^ + V\\£\\l2 + + |a|2 + |/3|2) . 

step 2 Computation of the localized virial identity. 



(4.23) 



'which is a consequence of the critical mass conservation 
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We now pick a large enough number ^ > and compute the locahzed virial 
identity: 



AVcf) 



Vee 



e + 



.v.(l). 



Ve 



y 



Vee 



(4.24) 



Now from ^J7} . IKm : 



1 /A 



A2 A V A 



+ b 



which yields: 



AVcp[^)-Ve?j=0{A\\e\\l,: 



We may thus rewrite (|4.24p as: 

y 



AV^ 



A 



• Vee 



j AVct) (I) • {dseiVe2 - a.esVsi)^ 
m A</. (I) {d,eie2 - a.eaei)) + 0{A\\e\\l,). 



We now inject the equation for e (|5.59p - (|5.60p and the estimates of Lemma l3.ll to 
derive: 



2A~ 



AV^ 



Vee 



AVcp 



A 



{M2{e)Ve2 + Mi{e)Vei) 



1 

2A 



(4.25) 

A0(|) {M2{e)e2 + Mi{e)ei) 



+ 0{A\\e\\]j,+V^ + V\\e\\L2 + \a\'^ + 



+ 



1 
A 
1 

2A 



1 



4^12 A 



j AV<P (V(Q2 + 2S2)e2 + 4V(Se)eie2 + V(Q2 + 2e2)e2) 



+ O {A\\e\\]j^ +V^+ V\\e\\L2 + jap + |/3|2) . 
step 3 Conclusion. 

Summing (I4.23P and ()4.25p and expanding the nonlinear term according to p.lSp , 
()3.16p . we obtain after a few algebraic manipulations: 



-Q,AQ 



2A 



AV(t) 



A 



Vee 



+ 



/vV(l),V.,W) + /|.P-ji,/A=*(l 



V^k{f)){a,a) 



3Q2e2 + g2^2 



{AVcp 



A 



y) ■ {3QVQel + QVQe^) 



+ OiAWeWl^,^. +V^+ V\\e\\L2 + \af + |/3|^ 



(4.26) 



41 



Now, the choice of the function (j) implies that: 

vV(f)(V..W)+/|.P-^/A^*(l)|.P 

> j e-'ii\Ve\'+ l\e\'+o(^M\l^. (4.27) 



Hence the quadratic form appearing in the RHS of (|4.26p is a smah deformation 
of the hnearized energy, and a standard locahzation argument^ together with the 
coercivity property (j3.30p . the orthogonahty conditions (|3.4p - (|3.8p and the conser- 
vation of mass (j3.23p ensure the existence of a universal constant c > such that 
for A > large enough: 



< l\e\'+ l\Ve\\~^-J{3Qhl + Qhl) ^^^^^^ 

+o(lkllHi(R^)+^')- 



Together with ()4.26p . ()4.27p and the bound (|3.17p . this yields the existence of con- 
stants c, c > such that for ^ > large enough and t{A) < t < close enough to 
0, we have: 



\ 4 



>|||a|2 + j + j \^je\^e'^^^+0{A\\e\\]j,+V^ + m% 

where we used the identity (jypQ, AQ) = —\yQ\^2- We may now fix A once and for 
all and (|iJ8]l follows. 
This concludes the proof of Proposition [ 



4.4. Convergence to of u in av^ay from the concentration point. We 

now slighlty change our point of view and consider the decomposition of u in original 
variables: 

u = Q + u 

with 

k{ait)y/^X{t) V m J (4.29) 
u(t x) = - e s e*7(t) 

Then the bound (|3.17p and the definition (|4.29p of u yield: 

II^(0IIl2(r2) < m and \\u{t)\\Hi^^2) < 1 for all to < t < 0. (4.30) 
Our aim is to improve the energy bound ()4.30p for the dispersive bound 

u{t) ^0 in as t ^ 0. 
The first step is dispersion away from the blow up point. 

Lemma 4.4 [H^ dispersion away from the concentration point). There holds: 

lim||n(t)||^i(|,|>3) =0. (4.31) 



19, 



see for example Appendix A in |20| 
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Proof of Lemma 14.41 

Let Q,u be given by (|4.29p . then u satisfies: 

idtu + Au = -R- k{x)u\u\'^ (4.32) 

with 

R = k{x) [(Q + u)\Q + n|2 - Q|Q|2 - ii|n|2j + id^Q + AQ + A:(x)Q|Qp. (4.33) 
step 1 Lfllf^^ bound away from the concentration point. 

We first claim as a consequence of the smoothing effect of the linear Schrodinger 
flow the space time bound: 

1-0 

/ \u{t)\'^ 3 dr < +00. (4.34) 

To wit, let the Fourier multiplier D = (1 — A)^/^. Let us start with treating the 
nonlinear term in (j4.32p . Let be the solution to 

idtC + AC = -k{x)u\u\'^ , C(0) = 0, (4.35) 

then from standard Strichartz estimates and smoothing for the linear Schrodinger 
flow, there holds: 

ll^l,.(,„,o),^i) + llCll^-(ao,o),H.) < ll^(^(-)-l-nil,|(,„,o),.f ) ^ 1 (4.36) 

where we used Sobolev embeddings and the L°°((to, 0), i?^) bound (jOOjl . 
We then let iv = D2 (^u — C) and x be a radial smooth function with 

^.r r for < r < 1, 

X{r) = / x"{p)dp and x"{r) =\ 1 for 2 < r < 4, 
•^0 [ ^ for r > 5, 

and such that 

Vr>0, ^<x"(r). (4.37) 



From (336]), (li^iD follows from: 

fO 

lL2(2<r<4) 
to y - - ' 



^M\h(2<r<4) ^ 1- (4-38) 



Indeed, from (|4.32p . w satisfies 

idtw + Aw = D^R, 

where R is given by (|4.33p . We compute the associated localized virial identity: 



(/ ■ ^**) } = - + ■ (4.39) 

= j x"{drwf + j jidrwf -^j AWw\'^ + K (^D^R, + Vx • Vw^ 

where drW = j.dgw. Let us now estimate the various terms in ()4.39p . The boundary 
term in time is bounded using ()4.30p . ()4.36p : 



<ll*li;,i <ll^^ll|^i(M2) + ||Cll|,i(M2)<L (4.40) 
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Similarily, 



< 



u 



m/2 



+ 



1^1/2 



< 1. 



(4.41) 



The R terms in (|4.39p are treated from Cauchy Schwarz using the space locaHzation 
of X- 



K ( D^R, ^w + Vx-Vw 



< 



\xD2R\ 



L2(r>l) 



Ax - 



-w 



+ 



L2(r>l) 



W 



L2 



< S I x"{drwf + \x\^\D\R\^ + \\w\\l. 



< 6 J x"{drwy + - J \xD-2R\' + l 

for some small 5 > 0, where we used (jilST]) and the bounds ()I30]l . (f06]l . 
We now claim 



(4.42) 



/ 

Jto 



\xD2R\\j^2 < +00 



(4.43) 



which together with (|i:40|) . (jOT]) . ([02]) injected into (|i:39|) yields (gSH]). 
Proof of (|4.43p : We use the usual symbolic calculus for the composition of pseudo- 
differential operators to obtain: 

D^\x\'^D^ = xDx + a^i{x,D) (4.44) 
where a_i(x, ^) is a classical symbol of order —1. In particular, a-i{x, D) is bounded 

4 

from L3 to L . Together with (j4.44p . this yields the bound: 

|2 I" Tih \ -ni Tj D\^ll7~iTT„Dll2 I ||DI|2 



\xDm\\l2 = {D^\x\^D^R,R) < \\D^xR\\l2 + 11^11^4. 



We change variables: 



R{t,x) 



A3 



S 



x-_aity\ 

m ' 



(4.45) 
(4.46) 



with 



A 



as (^^ XVk(a(t)) ^ \ 

- ^T-r^-+2TR^)r^-J+^^^- 

+ k{Xit)y + a{t)) [{Qv + e)\Qp + el"^ - QvlQvl'^ - e\e\'^] (4.47) 
which is well localized in y. Together with ()4.45p . this yields: 

\\xD"2R\\^2 < ^ [WD^Sh^ + \\d"2S\\l2 + 11511^4) . (4.48) 

We now explicitely expand the nonlinear terms in e in S and use standard commu- 
tator estimates together with the good localization in space of S{s, y), the bound on 
the geometrical parameters (j3.17p . (j3.20p . and the 0(\^) control in the construction 
of Q-p (|2:22D . ([2^ to conclude: 



||xD5i?||2, < ^ iVepe + llelli^ + + |a|4 ) < 1+ 
(OHIl now follows from (li:2T]l P1 



1 

A3 



iVere -/A + |k|| 



20recall that # = tt 

at 
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step 2 Strong convergence outside the blow up point. 

The strong convergence (|4.3ip now easily follows. Let a smooth cut off function 
■0 with -0 = 1 on |x| > 3 and -0 = on |x| < 2, then w = tpu satisfies the following 
equation: 

idtw + Aw = F — k{x)\ wfw, w{t) ^0 in as t 0, (4.49) 

with F given by: 

F = 2V^ • Vu + A'lpu - k{x)'ijj{l - V'^)|upM - tpR 

and R given by (|4.33p . (|4.3ip follows from: 

\\w\\L'=°i{to,o),m) ^ ^ as to^O. (4.50) 

Indeed, we write down Duhamel formula for ()4.49p and first derive a bound for Vw 
in L^((to, 0), -L^). Using the Strichartz estimates for the 2-dimensional Strichartz 
pair (4, 4) and the smoothing effect for the linear Schrodinger operator, we obtain: 

l|Vu;||L4((to,0),L4) + \\'^w\\L°°{{to,0),L^) (4-51) 

< ||V(fc(x)|«;pti;)||^4/3((i„,o)xR2) + \\\^\^Kmt,,o),Hhy 

The nonlinear term is estimated as follows: 

ll^(^(^)l^l'^)IL4/3((t,,o),Lf) ^ ll^^lli^((<o,0),L*)ll^lli4((to,o),L4) 
^ \NM\Li{{to,0),Li)\\u\\L°°((to,0),m)\\u\\L^((to,0),L^) ^ \\M\L°°{to,0)\\'^w\\L^i{to,0),Li) 

< ^\\^M\LH{to,o),L^) (4-52) 

for some small enough constant 6 > and where we used ()4.30p . The second term is 
estimated using the compact support property of Aip and (1 — ip'^) and ()4.34p . 

gas]): 

< o(l) as to^O. (4.53) 
(H3T]) . (|i32D and (|433]l yield: 

l|Vw||i4((tQ^o)xR2) + \\^w\\L°°{{to,o),L2) < o(l) as to 0, (4.54) 

and (j4.50p is proved. 

This concludes the proof of Lemma 14.41 

4.5. Convergence to in average of u in H^. Our aim is now to propagate 
the convergence of u away from the concentration point (|4.3ip to the blow up 
up region as well. Here the key is the refined bound (|4.2ip provided by the refined 
virial dispersive estimate (|4.18p and which implies : 



liminf^^=0. 

t^o X[t) 

The first step is to obtain a convergence in average in time. 

Lemma 4.5 [H^ dispersion in average in time). There holds: 

fO 
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Proof of Lemma 14.51 

Note that we may restrict the norm in (|4.55p to the region |x| < 3 in view of 
(lOTD . 



step 1 Morawetz identity. 

We first claim the fohowing virial type bound: 



r / \Vn\'<oi\t\) + \\umL^+ f 

Jt J\x\<2, Jt 



Il2 



+ / |Vepe ^ 



A2(r) 



dr. (4.56) 



Let X ^ smooth radial cut off function on ^ such that V^x is supported in |rE| < 4 
and is positive semidefinite in M^, and Vx(a^) = x on |x| < 3. We have: 

d 



2 / VMVn,W)-^ 



dt 



-2K 



(Au + A;(x)|i2ptt + R) 



Ax 

——u + Vx • Vu 
2 



1 



1 



Ax/c(a;)|n|* + o / Vx- V/c(x)|u 



L~,|4 



25R / i? 



Ax 



n + Vx • Vtt 



(4.57) 



where u, Q and R are defined as in (|4.29p (|4.33p . We estimate the various terms in 
the right-hand side of (|4.57p . Using (I4.30p . we have: 



1 
2 

< ll^ill? 



AxA;(x)|n|4 ) + - / Vx-VA;(x)|n| 



< 



u 



L2 



+ 11^1112 



IL2{R2) ^ ll"llL2(R2)ll"'llHi 
We now claim that we can treat the quadratic R type terms perturbatively thanks 
to the vanishing of the cut at order 1 at the origin. Indeed, we change variables 
according to (j4.46p with S given by (j4.47p well localized in y. We thus get from the 
bound on the geometrical parameters (|3.20p . Vx(x) ~ x near the origin and the 
O(A^) control in the construction of Q-p ^2^, ([2^ : 



(4.58) 



/ R 



Ax 



u + Vx • Vu 



1 

A2 



sft / S 



Ax(X(t)y + a(t)) 1^ ,,,, , ^ 
^ 'l —e + jVx{X{t)y + a{t)) ■ Ve 



< 1 



(4.59) 



We inject ()4.58p . (|4.59p into (|4.57p and integrate in time so that: 



\Vu{TfdT< 



< 



\x\<3 



9 / Vx • Vnn 



V2x(Vn,Vn) 



+ 



lap 



< 0(1*1) + l|f'(()llL» + 

where we used from ()4.30p : 



° lap 



A2(r) 

+ lk(r)|li.+/|V6(r)pe-^ 
A2(r) 



dT 



9 / Vx • Vuu 



< \\Vu\ 



l2|fIIl2 



< 



I^^IIl2, 



46 P. RAPHAEL AND J. SZEFTEL 

this is (|i36|) . 

step 2 Averaged in time dispersion. 

We now divide ()4.56p by \t\ and integrate in time. From the pointwise lower 
bound (|4.8|) and the dispersive bound (|I%TF I: 



1 /■0|a|2 + ||e(r)||2, +/|Ve(r)|2e-73 



\t\ Jt A2(r) 
Similarily, from Cauchy Schwarz, 



dr < 



o\a\^ + \\eir)\\l,+J\VeiT)\'e-VA 



A3(r) 



as t^O. 



\H'r)\\L2 



dr 



< 



° Mr)\\h 
A3(r) 



dr 



as t — ;> 0, 



and ()4.55p fohows. 

This concludes the proof of Lemma [53 



4.6. Control of the modulation parameters. We now claim that the dis- 
persion (|4.55p coupled with the conservation laws implies a sharp control of the 
modulation parameters and hence the derivation of the blow up speed. Similarly 
like in |8j, the key is to first derive convergence in a weak averaged in time sense, 
which thanks to the Lyapounov type dispersive control (|4.18p may then be turned 
into pointwise controls on the whole sequence in time. 



Proposition 4.6 (Pointwise dispersive bounds). There holds the pointwise bounds: 

+ \a\ + lkll//i(R2) < A^ 
b{s) 



X{s) 



1 



< y 



Kt) = -^ + 0{\t\^) for to<t<Q. 



(4.60) 
(4.61) 
(4.62) 



Moreover, there exists 70 G M such that: 

l{t)- 



^+lo + 0{\t\). 



(4.63) 



Remark 4.7. Note that ()4.60p implies in particular the zero momentum limit 
p.l2p . and ()4.60p improves the energy bound ()4.30p by a whole factor A. 

Proof of Proposition 14.61 

step 1 Control in average of a and (3. 
We claim: 



1 



lim — / ( — 

t-i>0 \t\ 



1 



A2 ^ A2 



d(T\dT = 0. 



(4.64) 



We start with a. Indeed, from (j4.2ip and the pointwise bound 
1 



\t\ 



° lap 



A2 



dT< 



° lap 



A3 



dr 



lap 



A 



dcr — )• as t 



2lRecall that 4f = 

at 



dr 
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which in particular imphes (j4.64p for a. 
We now consider /3. We have: 



a ■ (3 
IT 



/3 a a- /3 Xs 

2|/3|2 + (-6/? + co(a)A)-^ + 6^ 



(4.65) 



A 



a 



+ -i-2/3 -/S + (/?, + 6/3- co(a)A)--- -^ + 6 



A 



A. 



A 



a-/3 

IT 



= 2|/3|2+co(a)-a + 0(A3 + ||e||i2), 
where we have used ()3.20p in the last equality. Integrating ()4.65p between s and 



+00, we obtain: 

"+00 

|2, 



a{s) ■ (3{s) 
AC^T" 

< l«(^)l + 



r+oo /■+00 

/ co{a) ■ ada + / 0{X^ + \\e\\l2)da 

J s J s 



+ Ikllia) da 



X^{a)da 



where we have used ()3.17p which implies in particular that a(3/X — t- when s — )■ +oo. 
Dividing by |t| and using (I4.2ip yields: 



A2 



da < 



< 



< 



Hr)\ 
\t\ 

Hr)\ 
\t\ 

\t\ 



+ >^[l«(^)l' + lklli=^] +7/ ^(^)^^ 



da 

Ir W 
+ 0(1). 



[Ha)\' + \\e\\h] + o{l) 



Integrating once more and using Cauchy-Schwarz, this yields: 



A2 



2 \ 1 
da]<^ 



■dT+o(l) 



l«(r)P 
A3(r) 



dT +0(1) = 0(1) 



and (j4.64p follows for (3. 



step 2 Limit of on a subsequence as f — t- 0. 



and (ITO)) yield: 



1 /-o / 1 /-o / 1/3| 



lim — / I ■; — r 

t-5-0 



+ 



A2 ' A2 



+ 



A2 



da\dT = 0. 



In particular, this implies the existence of a sequence t„ — t- such that: 

A(t„) A(t„) A(f„) 
Injecting this into the conservation of the energy (|3.28p yields: 



lim 

n— >+oo 



Ktr^ 



We now observe from ()4.23p the bound: 



CI 



leP + lap 



+ 0(A2 



(4.66) 



(4.67) 
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and thus (|4.2ip ensures 



+00 



ds < +00 > 



Hence j has a hmit as s — >• +00, and from ()4.67p : 



Now from (fXTT)) . IKm . 



b_ 
A 



±— as t^T. 



A 



A 



+ b 



< A 



(4.65 



and hence the finite time blow up assumption together with A(t) > imply: 

lim ^ = > 0. (4.69) 

Injecting this into (j4.68p and integrating in time yields in particular the pseudo 
conformal speed: 



A(t) = 11(1 + 0(1)) as t^O. 
Co 



(4.70) 



step 3 Improved bounds. 



We now claim that the knowledge of the asymptotic limit (|4.69p coupled with the 
monotonicity (|4.18p and the conservation of energy implies a spectacular improve- 
ment on the bounds of ||e||j:^i. 

Indeed, we integrate the local virial identity (|4.18p between t and t„ and let t„ — )• 0. 
The boundary term in t„ is estimated using (j3.17p . (j4.69p : 



bjtn) 

A(t„) 



4 ^ + 2A^ 



AV(t) Ve(t„)e(t, 



1 

Co 



(4.71) 



as tn — 0, and we thus get: Vs > 0, 



(f + Mh} + 0{V' + + \\e\\l,)) da. 



X{s) Co, 

+00 



which together with (|3.17p implies: 

/+00 -j^ 



< 



w- ) + Mm 



X{s) Co 
We now recall (|3.29p which implies: 



+ 



(A3 + 



+ k 



H 



)da. 



(4.72) 



|2 + lap + ||e||^i 



Dividing by A^, we obtain: 



A2 ^ A2 ^ 



|2 

\m 



A2 



A2 



+ A^. 



(4.73) 



(4.74) 
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Finally, multiplying (|4.72p by + adding to (|4.74p . and noticing that the 
terms — ^ cancel each other yields: 



A2 



/• + 00 /■ + 00 



< + / (A^^ + + \\e\\jj.)da < + X'da <\t\<X 

J S J S 

where we used (|3.17p . (|4.70p . This yields the improved pointwise bound 

m' + \\e\\h < 

which reinjected into (|4.75p yields: 



(4.75) 



/ + 00 -j^ 
j{\a\^ + Mh}da 

/•+00 r+oo 

< A2 + / (A3 + |/3|2 + \\e\\jj,)da < X^ + X^da 

J S J S 



— h - — ■ h 

A2 ^ A2 ^ 

f + OO 



< y+ xdT < X' 



(4.76) 



where we used (j4.70p again. This concludes the proof of the improved bound (j4.60p . 
We now integrate the localized virial identity (j4.18p between s and +oo and estimate 
using Km . K76\i : 



b{s) 



1 



X{s) Co 



/H-oo -j^ r-\-oo 
J W + MI2} da + ||e||Hi(R2) + J X^da 



< A^ 



and KUB is proved. Next, from (H^QD, (OOjl . dMIJ: 

i^ + b = 0{X^) and thus {X)t + ^ = 0{\t\^) 
A Co 

which yields ()4.62p by integration in time. This eventually implies using also p.20p . 
(TO : 



(4.77) 



d 



r^2 



1 

A2 



Is 



clx^ 



t2 



1 

A2 



2x2 



t2 



< 1 



and Km follows. 

This concludes the proof of Proposition 14.6 



Remark 4.8. In view of ()4.60p and (|4.6ip . we have: 

In particular, the approximate solution Q-p is of order five -see (j2.9p - as announced 
in Remark 12.21 Also, in view of (|4.60p . we have: 



(4.79) 



which together with (|3.20p implies: 



a + 



+ \Modit)\<X\ 



(4.80) 
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5. Uniqueness 

We now have obtained the exact blow up speed of critical mass blow up solutions 
which corresponds to the crucial relation 6 ~ A and the dispersive behavior 

n — > in as t — > 0. 

Let now Uc be the critical mass blow up solution given by Proposition 13.61 which 
blows up at T = and a* = 0, with energy Eq and phase parameter 70 given by 
(j4.63p . We need to prove that 

u = Uc- 

The proof proceeds in two steps. We will first show that the refined estimates of 
Proposition 14.61 together with the backwards propagation of smallness of Lemma 
13.51 implv the strong convergence 

lie — u — )• in as t — t- 0. 

We then show that the a priori estimates obtained for u — Uc are strong enough to 
treat perturbatively the growth induced by the null space of L when linearizing the 
equation close to Uc and running the energy estimates of Lemma 13.31 Both these 
steps require having sufficient a priori decay estimates and in particular require the 
construction of an the approximate solution to at least the order O(A^) for the error. 

5.1. convergence to the critical element. We claim the following dispersive 
property which somehow collects all previous estimates on the solution and is the 
key to the proof of uniqueness. 

Lemma 5.1 (H^ convergence to the critical element). There holds the strong con- 
vergence at blow up time: 

u — in H^. (5.1) 

More precisely, 

||V(7x-n,)||i2 + ^i^|^^<|t|3 as t^O. (5.2) 
Proof of Lemma 15.11 



step 1 Backwards propagation of smallness and improved bounds on the solution. 

We first claim that the bounds of Proposition 14.61 coupled with Lemma 13.51 imply 
a refined bouncF^: 

\\m\L^<^'{t), l|Vn(t)k2 <A3(t). (5.3) 
Indeed, we decompose u according to the geometrical decomposition (j3.ip . We then 
let an increasing sequence of times such that i„ — )• 0. Observe from the conservation 
of mass and Proposition 14.61 that the assumptions (|3.53p - (|3.55p of Lemma 13.51 are 
satisfied at any t„. In particular, there is a time to < such that we have ()3.56p for 
any to <t < tn- 

iNmwh + < iivn(t„)iii. + + m 

We now let n -)> +00 so that ^^m\i yields (fO|) . 

step 2 Comparison between the modulation parameters of u and Uc- 



'Note that this bound is the same like for Uc, see (|3.82p 
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Let b, A, a, (3, 7 denote the modulation parameters of u, and Ac, Oc, /3c; 7c denote 
the modulation parameters of Uc- We claim: 

|7 - 7c| + ^yp^ + + |fe - bcl + 1/3 - /3c| < (5.4) 

The proof of (|5.4p is a consequence of the modulation equation and the improved a 
priori bound (|5.3p . and it is postponed to Appendix C. 



step 3 Comparison between u and Uc- 

(j5.2p is now a simple consequence of (|5.3p . (j5.4p . Indeed, 

^^p^ + ||V(n-nc)||L. < ^ + I|Vu||^. + ^|^ + ||V^c||l^ 

+ ^p|^ + I|v(Q-Qc)IIl- 

Now (|5.4p yields after a simple computation: 

< |^|^ I|V(Q-4)IIL2<|^|^ 

and (1531) now yield CT . 

This concludes the proof of Lemma 15.11 

5.2. Energy estimates for the flow near Uc- Let us now decompose: 

Here we do not impose modulation theory and orthogonality conditions on e. In- 
deed, one should have in mind that we do not have uniform well localized bounds 
on the Uc part of Uc and that the dtUc part is in particular controlled only through 

3 

the H2 bound (|3.82p . We however claim that the a priori estimate from ()3.82p . 

II^IIl2 < AMI^IIhi < A3 (5.6) 

is enough to treat the instability generated by the null space of L perturbatively. 
Let 

Nit) := sup ( ||^(r)|||. + ) , (5.7) 

t<T<o \ A^(r) J 

and 

Scal(t) = {ei,Qf + {e2,AQf + (ei,\yfQf + {e2,pf (5.8) 
+ (ei,yQ)' + (e2,VQ)2. 
We first claim the following energy bound: 
Lemma 5.2. There holds for t close enough to 0: 

N(t)<sup^+["^4r. (5,9) 

t<r<o A-(r) Jt A^(r) 

Proof of Lemma 15.21 

It is a consequence of the energy estimate (|3.38p together with the a priori bound 
(13:82]) . 
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step 1 Application of Lemma 13.31 
Let I{u) be given by p.37p . we claim that: 



A3 V \l ) ^ dt' 

Indeed, let us apply Lemma [3.31 with w = Uc = {uc)i +i('Uc)2) then the bound (|3.32p 
holds from (j3.82p and ip given by (|3.3ip is identically zero. Hence (|3.38p becomes: 

^ = (^j k{x)ulu ^ - K k{x)dtUc{2\u\'^Uc + v?u^) 

+ T2 W / ^ <^ — (Vn,Vn 



+ o(^^ + Mjj?j. (5.11) 
We consider the first two terms in the right-hand side of (jS.lip and expand Uc = 

Qc + 

-^9 (y Hx)ult^^ ~ ^ (/ ^(^)^t^c(2^S^n~fl2^^ (5.12) 



- K / A;(x)atnc(2|nPQc + ii^Qc) - / k{x)dtUci2\u\'^Uc + u'^Uc) 



Arguing like for the proof of (|3.72p . we may rewrite the first two terms in the 
right-hand side of (|5.12p as: 

-^9 (/ ^(^)^'^') " ^ (/ Hx)dtQc{2\ii\^Qc + (5.13) 



Keep in mind that we do not have satisfactory well localized bounds in Uc, and the corre- 
sponding terms will be treated using the smallness (|3.82|l 
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where Qc = T,c + iQc- For the next two terms in the right-hand side of (|5.12p . we 
use Sobolev embeddings and (|3.82p to obtain: 



^ l^llQc||L'^{M2)ll^f. 



'c||L2(R2)||'""^ 



+ ~ ll^cll j;,4 



+ l|f^iQc||L°°(R2)ll^^c||L2(iR2)||tt||^4(K2) < 

where we used the bouncrj 



xr 



U\ 



L4(M2) 



c||L°° ^ 



A2 



+ u 



(5.14) 



The last two terms in the right-hand side of (|5.12p require using the equation sat- 
isfied by Uc- 

idtUc = -Auc - k{x){\uc\'^Uc - IQd'^Qc) - i^c (5.15) 
where if^c is defined by: 

i^c = idtQc + AQc + k{x)\Qc\^Qc. 
Expanding this term hke for ([STTi]) . we have using l^320h . (|332]), ^Ml, 

Modcit) + \ac\^ + Xl , ^2 



A? 



< A^ 



Using this together with (j5.15p . integration by parts, Sobolev embeddings and the 
H2 bound (|3.82p now yields: 



/ k{x)dtUci2\u\^Qc + u'^Qc)j -^yj Hx)dtUci'^\u\'^Uc + u'^Uc) 

2„. |2 



\\Uc\\ 3 



\\2\u\^Qc + v?Qc\\ „i + ||2|nruc + w Uc|| „i 



+ \\k{x){\Uc\^Uc - \Qc\^Qc) + V'c||l2 IIQc||L°o(R2)||u||i4 + || || L6 ll^^llie 



< 



A? 



+ \\u\ 



(5.16) 



We now consider the last term in the right-hand side of ()5.1ip and compute: 
= 5ft ( /" AV(I) 



V AX, 

X — a 



A;(x)(2|?2| Uc + u Uc)Vuc 



(5.17) 



AX, 



k{x){2\u\^Qc + u^Qc)VQc Error 



24- 



The worst term is generated by the phase |(7c)t| ^ 



54 



P. RAPHAEL AND J. SZEFTEL 



with from Sobolev embeddings and (j3.82p : 

' X — a, 



I Error I 



+ 



+ 



3? \ AV^ 



y ASJ(\> 
X — ac 

X — OLc 



A\r 



k{x){2\ii\'^Qc + u^Qc)Vu, 



?.2, 



AXr. 



< (||VQc||l°°||'Uc||l2 + IIQc||L°°||Vnc||L2)ll'" 



?,l|2 
L4 



+ ll^'c||L6||Vnc||i2||n||^6 
Collecting the estimated (|5.1ip - (|5.18p yields: 



(5.18) 



dX be 



A2 



X — Up 



+ ^{ I V^</>(-^)(Vn,Vn) 



4^2 



+ o 



(/(av*( 



X — Of 

^4A7 



X — 



X — a, 



k{x){2\u\^Qe + rQc)-VQ, 



li=,l|2 



A2 



We now use the uniform proximity of Q-p^ to Q and the coercitivity property (|3.30p 
to conclude like for the proof of (|3.72p : 



6c 



\y\ 



iVe^e ^ + / |er + 0(Scal(t)) 



+ 



n 



II ~ ||2 



which implies (jS.lOp . 
step 2 Coercivity of X. 
We now recall from (j3.37p the formula: 



At) = ^ I + ^ I I k{x)\uc + u\'' + ^ I k{x)\u. 



1 f\W 1 



A2 



1 



1 b. 



+ I kix)\ucr{uc)iui+ I k{x)\Ue\^{Uc)2U2 + -Y^\ I 



X — Oc 
AXr 



■ Vuu 



Expanding Uc = Qc + Uc and arguing like for the proof of (jS.lOp . we get using (|5.6 
the rough upper bound: 



1X1 < 



_ Mt)\\m + ;^2(^) ^0 as t^O. 



(5.19) 



We now claim the lower bound: 

1 



X(t) > — [(L+ei,ei) + (L_e2,£2) + o(lkll^i)] 

Scal(i) . 



> 



2A2 

c 

a! 



|p||2 



(5.20) 
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The proof is very similar to the one of (|3.65p using also the control of interaction 
terms similar to (|5.14p . (|5.16p . (jS.lSp . The details are left to the reader. 
Integrating ()5.10p from t to using the boundary condition ()5.19p and the lower 
bound (fOOjl now yields (fOjl . 
This concludes the proof of Lemma 15.21 



5.3. Control of the scalar products and proof of Theorem 11.31 It now re- 
mains to control the possible growth of the scalar product terms in (|5.9p . We claim: 



Lemma 5.3 (A priori control of the null space). There holds for t close enough to 
0: 

Scal{t) <\t\^t\'^N{t). (5.21) 
Let us assume Lemma 15.31 and conclude the proof of Theorem 11.31 

Proof of Theorem 11.31 

From (|5.9p . (j5.2ip and the law Ac ~ \t\, we have for t close enough to 0: 
N{t) < \t\'^Nit) + J' ^dr < \t\'2N{t) 

and hence N{t) = for t small enough. From the definition (|5.7p of N, this yields 
u = Uc and concludes the proof of Theorem 11.31 

Proof of Lemma 15.31 

The proof follows by deriving the null space close to Q-p^ to sufficiently high 
order and reintegrating the corresponding modulation equations from blow up time. 
The worst behavior is on the even terms where the modulation equations are a 
deformation of the ones for L+, L_, and roughly correspond to the system of ODE's: 

(e2,AQ). = 2(ei,Q), (ei,|y|2Q), = -4(e2,AQ), {e2, p)s = -{ei,\y\^Q) (5.22) 

with initial degeneracy provided by the norm conservation law and the a priori 
bound (ESD: 



\{euQ)\< / \e\'<\\e\\L^X{t)^/m<XHt)mVm- 



The control of the worst paramater (related to the phase) requires 

p+oo r+oo r+oo 

\{^2,p)is)\< dsi dS2 \i£l,Q){s3)\dS3 

J S J Si J S2 



< x{s)y^) p ds, ds, xtds3 < Mx^. 

J S J Si J S2 ^ 

This implies (|5.2ip and explains why we needed a small enough estimate ||e||j;,2 ^ 



in (I521PI, (ESI). 



Let us now implement the above strategy which requires being careful with respect 
to polynomial losses, and in particular demands a high order approximation of the 
null space close to Qp . 

step 1 Approximate equation in conformal variables to the order O(A^). 



25 



In fact, any bound ||e||j^2 < Xi^^ for some r; > is enough 
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Let v,w he defined by: 



1 1 / X — ac{t)\ (f) , , , , i hsMl -i B -ij 



uit,x) = — — -— TYTTT^' t, e'^'^^ ^ w{s,y) = v{s,y)e 

(5.23) 

then w satisfies tlie equation: 

{{Pc)s + bcl3c) ■ y + lKl3c ■ -J-, 77^ \ W + W\w\ 

k{ac{t)) J k{ac{t)) 
+ b,) (aw + (-ibj-^ + if^cv) w] + ((7c). - m^)w 



Ac 7 V V 4 

Starting with Uc, we also define Vc and Wc{s,y) = Vc{s,y)e^ 4 *Pcy_ -yy^ jg^ 
u = Uc + u and define: 

V = Vc + e, w = Wc + e, i.e. e = ee* 4 ^P'^'^. (5-24) 
Since satisfies (jl.ip . e satisfies: 

+ - e + ((6,), + 62)^ _ /((/3^)^ + b,p,) . y + iAc/3c • ^^^^^ 

4 [ k{ac{t)) 

k{Xc{t)y + ac{t)) 2 I |2n N |«|2n~ 

('^)^+6c^ ^e+f-.6c^ + ^/3cyV~ 



gJOD-gSH), Ml and (fOSll yield: 

+ A.- - e - (co(ac)Ac + /JsA^) • + M^4^i±^(^|^|2 _ ^^|^^|2) 

A;(ac(t)) 

= 0{Xtil + \y\')e + Xtil + \y\)Ve). (5.26) 

Let e = £1 + ie2 and Wc = {wc)i + i{wc)2- We now expand the nonlinear term in 
()5.26p as well as Wc = P-Pc + and the expansion of the approximate solution (|2.7p 
to derive the equation at order O(A^): 

- idsi + M^^\e) + (co(ac)Ac + P^X^) ■ ye = -t/; (5.27) 

where M^'^^ is the expansion of M at order 4: 

M(4)(e) = )(ei) + iM^ ^(ea), (5.28) 

Mi^\e) = -Asi + £1 - TsQ' + eQTa + ^V^k{0){y,y)XlQ^ + 3V^k{0){y, ac)XcQ^ 



+QQT3 + ^V^k{0){y,y,y)XlQ^ )ei - 2QS-ie2, 



(5.29) 
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M^''\i) = -Ai2 + 82 - {q^ + 2QT2 + ^V2fc(0)(y,y)A2Q2 + v2A:(0)(y, q,)A,Q2 

(5.30) 

and where the remainder tp satisfies: 

iP = (A^(l + |yp)e + A^(l + |y|)Ve + e^e + ele + u;^^^ + e^) . (5.31) 

step 2 Approximate nuh space. 

Let /(s, y) = 0(e~'^l^l) be a smooth weh locahzed slowly time dependent function, 
then yields: 

^ {9(eJ)} = M(4)(/) - + (co(ae)Ae + hXl) ■ y/) + 0((V^, /)) (5.32) 

with 

KV';/)! ^ Ac||e||L2(]i|2) + ||ec||L2(R2)l|e||L2(R2) + ||£'c|li4(R2) ||£||l2(]r2) 
+ ll'?c||L4(IR2)||e||i4(K2)||e||i2(]R2) + ||e||i4(]R2)||e||L2(R2) 

^ ^^£11^1^2(^2). (5.33) 



'c I 

We now claim that we can find some real valued smooth well localized functions 
M (y) ,B2{y),D3{s,y),E2{s,y),Es{s,y),F2{s, y) with 

A2.B2 = 0{\le~'\y\), E2,F2 = 0{Xle~'\y\), i?3, ^3 = 0(A^e-*l) 

such that we have the following approximate null space relations: 

mW(VQ + XIA2) - ids{VQ + XIA2) = -aiXlyQ + ©(A^e"!^!), (5.34) 

M^^\i{yQ + XlB2))-ids [i{yQ + XlB2))\ = -2i(VQ + A^^a) + 0(A?e-l^l), (5.35) 
for some universal constant 

ai > 0, (5.36) 

and for the even part: 

M^^\APv^ + E2 + + iD-i) - ids [^Pv, +E2 + E3 + iDa] 
+ (co(ac)Ac + P3XI) ■ y [APr, + E2 + E3 + iD^] 

= -2Pp^ + Xca2{ac)-yQ + X^,a3-yQ + 0{Xie-\y\), (5.37) 
for a linear map 02 on M? and a vector 03 in M?, 

M(4)(i(|y|2p^^ +F2)) - ids [i{\y\'Pv.+F2)] 

= -4iiAPv, + E2) + OiXle'\y\), (5.38) 

and 

M^^\p) = \y\^Q + 0{Xle-\y\). (5.39) 

The algebraic derivation of (Oil) . (lOSD . (lOeD . (lOTD . (lOHD . dOOD is done in 
Appendix D. 

step 3 Control of (ei,yQ). 
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From (lOSD . (lOD . (fOSjl . let A = VQ + A^Aa, B = yQ + A^^a and rewrite 
([5:32]) with / = ^ and then f = iB to get: 

(ei,S), = -2(62,^) + 0(A3||e||i2(j 
(£2, = aiA2(ei, yQ) + 0(A3||e||i2( 

Now, using the fact that ^ = ^ together with (|4.62|) yields: 



(5.40) 



A.(.) = ^ + 0(^). (5.41) 



([530]), ([OT]) and the fact that B = yQ + 0{\l) imply: 

(ei,S), = -2(e2,A)+Fi, / 



We conclude from Lemma 15.41 with = aiC^ > from (I5.36P that 

r+°° Ma) 

\{£i,B)\ + \{e2,A)\< / ^2^||e||i2(iR2)(ia. 
Using the fact that A = VQ + 0{s~^) and B = yQ + 0{s-^) yields: 

|(ei,2/Q)| + 1(62, VQ)| < ^ 
This yields together with (fOill . (fSlSHIl . (fSlMl) and the definition (fSTll : 

\{£i,yQm\ + \{e2,VQm\ < |i|||^||i2(K2) + ^ ln(|cT|)||^||i2(R2)da 

< |t|^|t|v^ 
and (j5.2ip is proved for the odd directions. 

step 4 Control of {ei,Q), (e2,AQ), (ei, |ypQ), (£3, p). 

We first have from the conservation law, the critical mass assumption and 



Hence from (13:821) : 



^{£,PvJ = 2?li{£,w^) + 0{Xt\\e\\L2) = 0{Xi\t\^/N{t)). (5.42) 
Using (lOHD . (lOTD and (l5:32]l with / = AP-p^ + E2 + E3 + iL>3, we have: 

= 2K(e,P-pJ - Aca2(ac) • (ei,y(5) - A^ag • (ei,yQ) + 0(A^||e||i2(R2)). 
(lOTIl and (fOajl yield: 



|a(e,APp, +£;2 + ^3 + ^^3)1 
Together with ([02]) . (fSliS]) and the definition (fSTll of iV(t), we obtain: 
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Using (lOSD . (lOHD and (lOSjl with / = i\y\'^Pv, + ^-^2, we have: 

- (5ft(e,|y|2Pp^ + F2))^ = 49(e,APp, + i?2) + 0(A3||e||i2(K2)). 
([OT]) and (lOSjl yield: 

3f?(e,|y|2Pp^+F2)| < (^|!3(e,APp, + i?2)| + "^"^3^'^'^ ) ds. 

Together with ([5li¥|) and the definition (jST]) of iV(t), we obtain: 



5R(e,|y|2Pp^+F2) < |t|5|t|0V(t). 
Using (lOHIl . (ICTH and (f02]l with / = p, we have: 

(e2,p), = -(ei,|?/|^(5) + 0(A^||e||i2(B;2)). 
(|5:iT]) and (fOHjl yield: 



l(e~2,p)| 



l"^{\ier,\y\'Q)\+^^^)ds. 



Together with (f5^ . the fact that P-p^ = Q + 0(t2) and the definition §. 
we obtain: 

\{e2,p)\<\t\'^\t\Vm. 

Finally, (f05|) . dEHD, (IS3ZD, ([^30]) together with the fact that P-p^ = 
imply: 

+ 1(62, AQ)| + 1(^1, |y|'g)l + 1(^2, P)| < \t\^\t\y^ 
and ()5.2ip is proved for the even directions. 



(5.45) 
(5.46) 

(5.47) 
(5.48) 

(5.49) 

of iV(t), 

(5.50) 
Q + 0{t^) 



Appendix A: Computation of the modulation parameters 

This Appendix is devoted to the algebraic computation of the modulation param- 
eters in the setting of a decomposition (j3.ip satisfying the orthogonality conditions 
()3.4p - ()3.8p . We in particular aim at deriving the estimates ()3.20p . ()3.2ip . Similar 
computations can be found in [29], [30]. We will in fact prove a slightly more general 
result which will be used in Appendix C. 

Let Mod{t) be the vector of modulation equations given by ()3.19p . we claim that 
we can find a polynomial vector IZ with 

\R{n < V{\a\' + |/3|2) + |P|4 + (ft - ^) irl' (5.51) 

such that: 

Mod{t) = n{V) (5.52) 
+ 0{V' + [V' + \Mod{t)\]\\e\\L2 + \\e\\l2 + \\e\\%^ + \\\u\\l2 - WQh^l) 
and this implies in particular (|3.20p . 



step 1 Inner products. 
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We compute the inner products needed to compute the law of the parameters 
from the Qv equation ^(TTI^ where Mi,M2 are given by (j3J3]) . ^J^: 

- i-M2{e) + 6Aei - 2/3 • Vei, VG) + (Mi(e) + bAe2 - 2/3 • Ves, VS) 

{-M2{e) + bAei - 2/3 • Vei,yS) + (Mi(e) + bAe2 - 2/3 • Ve2,ye) 
= 0{V^e\\L2), (5.54) 

- (-M2(e) + bAei - 2/3 • Vei, AO) + (Mi(e) + bAe2 - 2/3 • Ve2, AS) 

= -2^ie,Q^) + 0{V^\\e\\L2), (5.55) 

(_M2(e) + 6Aei - 2/3 • V^i, + {Mi{e) + 6Ae2 - 2/3 • Ve2, lypB) 

= 0(P'|k||L2), (5.56) 

- (-M2(e) + 6Aei - 2/3 • Vei,p2) + (Mi(e) + 6Ae2 - 2/3 • Ve2,Pi) 
= 0{V^e\\L2). (5.57) 
To prove (|5.53p - (|5.57p . remark using the definition (|2.2ip of Q-p that: 

= Qe-^^^+*''-^ + 0(p2). (5.58) 

|^|2 

Thus, it suffices to prove (f533]l - (15371) with Q-p replaced by Qe-'''—'^''^-y. This is 
done in |30j for (|5.53p - ()5.56p . Let us briefly sketch the proof for ()5.57p . From: 

Mi{e) = -Aei + ei - {{Q^ + 2^'^)ei + 2See2) + 0{Vh), 

M2{e) = -Ae2 + £2 - {{Q^ + 2e^)e2 + 2Seei) + 0{Vh), 

we have: 

- {-M2ie) + bAei - 2/3 • Vei, /Oa) + (A/i(e) + bAe2 - 2/3 • Ve2, Pi) 
= {ei ,-Api+pi- (Q2 + 2S2)pi - 2m p2 + bAp2 - 2/3 • Vp2) 

+ (£2, -Ap2 + P2- {Q^ + 2e^)p2 - 2Se/9i - bApi + 2/3 • Vpi) + 0{rh). 

Next: 



S = Q + 0{r'), G = (^-6^ + /3 . yj Q + o{r'), 
Pi=P + 0{V^), p2 = i-b^^ + /3 • y V + 0{V^)^ 



from which: 

- {-M2{e) + 6Aei - 2/3 • Vei,/)2) + (Mi(e) + 6Ae2 - 2/3 • Ve2,pi) 
= (ei,L+(p))+ |^e2,(-6^ + /3y)L-(p)-2(-6| + /3) • Vp - /5(-6) 

-2 (^-6^ + /3y) - hAp + 2/3 • Vp^ + 0{V^\\e\\L2) 



= {ei,L+p) + \e2,\-b'-^ + (3-y'jL+p'j+0{V'\\e\\L2) 

= (ei,|y|2s) + (e2,|ype)+0(p2||e||^,). 
The orthogonality condition (j3.7p now yields (|5.57p . 
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step 2 Simplification of the equations. 

To compute the modulation equations driving the geometrical parameters, we 
first simplify the equations (|3.1ip and (|3.12p using the explicit construction of Q-p 
as follows: 

(bs + b^)dh^ + i/3s + b(3- co(a)A - /JgA^ - 13^X^)8^^ 
+dsei - M2{e) + bAei - 2/3Vei 

h + 6^ AS + {% - |/3|2)e + (y - 2/3) VS + - R2{e), 



(5.59) 



and 



{bs + b^)dbQ + (/3, + 6/3 - co(a)A - /^sA^ - /34A^)a^e 
+dse2 + Mi(e) + 6Ae2 - 2/3Ve2 



y + 6)Ae-(7. 



|2)S+(^-2/3)VS-5R(V^) + i?i(e) 



(5.60) 



where the remainder ijj has he following form from (|2.23p . ()2.9p : there exists Fp{y) 
and a matrix S'p{y) polynomial in V and smooth and well localized in y such that 

i;{y) = Fr{y)+Sv{y)Mod{t)+0 {V^ + Mod{t)) [\e{y)\ + |(1 + |y|)Ve(y)|]+0(P^)e-^IJ'l 
with 

|5,'5p(y)| <|P|e-^'=|j^l, (5.61) 



3 \ g-Cfc|j/| 



(5.62) 



step 3 The law of b. 



We take the inner product of the equation (j5.59p of ei with —AQ and we sum it 
with the inner product of equation (|5.60p of £2 with AS. We obtain after integrating 
by parts: 

|2 



(6, + 6^)<j (^Q,AQ 



2K(e,Qp) 



= (i?i(e),AS) + {R2{e), A@) +ni{V) + Si{V)Mod{t) 
+ OiV"^ + \Mod{t)\)\\e\\L2+0(V^) 

with TZi polynomial satisfying (|5.5ip . Si polynomial in V of order 1, and where 
we have used the orthogonality condition (|3.6p and the computation of the inner 
product (|5.55p . We now rewrite the conservation of mass (|3.22p : 



2»(e,Qp) 



\e\' + 



+ k{a 



\u\ 



\Qj,\'-k{a) / Q 



+ n2iV) (5.63) 



where from ()2.25p TZ2{V) is polynomial and satisfies (|5.5ip . This yields using the 
definition (l3T8]) : 



{bs + b')\ ( ^Q,AQ 



V^k{0){a,a) 



WyQW 



L2 



{bs + b^ - (io(a,a)) 



= - j \e\'^ +n^{V)+Si{V)Mod{t) + (i?i(e),AS) + (i?2(e),Ae) 
+ O {{V^ + \Mod{t)\)\\e\\L2 + \\\u\\l2 - WQWlA+V') 



(5.64) 
(5.65) 
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with 7^3 polynomial satisfying (|5.5ip . 
step 4 The law of A. 

We take the inner product of the equation (|5.59p of ei with and we sum 

it with the inner product of equation (|5.60p of £2 with jypG. We obtain after 
integrating by parts: 

= (^^ + {{AQ,\y\'^Q)} +TZ4{r) + S2{r)Mod{t) (5.66) 

+ O ((p2 + \Modit)\)\\e\\L2 +V^ + \\e\\l2 + ||e||f^i) 
with 7^4 polynomial satisfying from (j5.62p : 

17^4(^)1 < r{\a\^ + |/3|2) + + (& - ^) 1^1', (5.67) 

^2 polynomial in V of order 1, and where we have used the orthogonality condition 
()3.7p and the computation of the inner product ()5.56p . 

step 5 The law of a. 

We take the inner product of the equation ()5.59p of £1 with yT, and we sum it 
with the inner product of equation (|5.6Up of £2 with y@. We obtain after integrating 
by parts: 

= (^-2p){{yQ,VQ)}+n5{V)+S3{V)Mod{t) (5.68) 



A 

+ O {{V^ + \Mod{t)\)\\e\\L2 +V'^ + \\e\\l2 + \\e\\%i) 

with 7^5 polynomial satisfying (j5.5ip . ^3 polynomial in V of order 1, and where 
we have used the orthogonality condition (j3.5p and the computation of the inner 
product (|5.54p . 

step 6 The law of /?. 

We take the inner product of the equation ()5.59p of £1 with — VB and we sum it 
with the inner product of equation (|5.60p of £2 with VS. We obtain after integrating 
by parts: 

-iPs + bf5- co{a)X - /33A3) {{yQ,VQ)} = TZ^^iV) + +Si{V)Mod{t) 
+ 0{{r^ + \Mod{t)\)\\e\\L2+r'^ + \\e\\l2 + \\e\\%^) (5.69) 

with TZq polynomial satisfying ()5.5ip . ^4 polynomial in V of order 1, and where 
we have used the orthogonality condition ()3.4p and the computation of the inner 
product (|5.53p . 

step 7 The law of 7. 

We take the inner product of the equation (|5.60p of £2 with p. We obtain after 
integrating by parts: 

-{bs + b^) I f ^^Q, p]] = {% - |/3|') {{Q, p)} + TZriV) + S,{V)Modit) 



_ 4 

+ 0{iV'^ + \Mod{t)\)\\e\\L2+V'^ + \\e\\l2 + \\efHi) (5.70) 
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with TZj polynomial satisfying (j5.5ip . polynomial in V of order 1, and where 
we have used the orthogonality condition (j3.8p and the computation of the inner 
product (|5.57p . We now inject ()5.64p and the definition (|3.18p and get: 

(7, - \f3f + di{a,a)){{Q,p)} = lZs{V) + SQ{V)Mod{t) 

+ o + \Modm\\4L^+'P^ + M\h + + \M\l^ - WQWlA) 

with 7^8 polynomial satisfying (|5.5ip . Sq polynomial in V of order 1. 
step 8 Conclusion. 

Putting together (|5.64p , (j5.66p , (I5.68P , (j5.69p , (j5.70p and estimating the nonlinear 
interaction terms using Sobolev embeddings yields: 

{A + Sj{V)} Mod{t) = ng{V) 

+ O {[V^ + \Mod{t)\]\\eh2 + Ml, + \\ef^, + \\\uh2 - WQh^l) 

for some polynomial vector TZg satisfying ()5.5ip . some 0(1) invertible matrix A and 
some polynomial matrix Sj(V) of order 1 in V. Inverting this relation to compute 
Mod{t) and computing the Taylor expansion of (^ + 57)"^ to sufficiently high order 
yields ([532]) which implies (f3l^ . 

We now now inject (|3.20p . (|5.67p and the almost critical mass assumption (|3.3p into 
(|53HD and (fOT]) follows. 

Appendix B: A priori bounds on an ODE 

In this appendix, we study an ordinary differential system which will intervene 
in Appendix C and Appendix D. 

Lemma 5.4. Let Z = (Zi, Z2) : M — t- satisfying the following ordinary differen- 
tial system: 

/ -2 \ 

Zs=\ ? n ]Z + F (5.71) 



s 



2 



where ? > is a constant and where F = [Fi, F2) : M — ?• M^. Assume also that 

lim Z{s) = and \F{s)\ < \. (5.72) 

s— >-+oo s 

Then, we have the following estimate for Z for s > 2: 

/+00 
{\Fi{a)\ + \F2{a)\a)logada. (5.73) 

Proof of Lemma 15.41 We start with the homogeneous differential equation. We 
look for Zfi solution of: 

{Zh)s = I ~^ ^ Zh. (5.74) 

We have: 

where : M — t- satisfies the following second order differential equation: 

. 2? 
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We look for solutions of the form z = s"^ , which yields the equation: 

t?2 - + 2? = 0. 
We find a basis , z_ with wronskian 



W{z+,z-) = -^z+{z^)s + ^iz+)sz- 



as follows: 



if < ? < ^, then i?± = ^ ^ z± = /± and W = VT^, 

8 2 

if ? = -, then , 2;+ = ^/slogs, Z- = y/s and W = 
8 2 



and 



if ? > ^, then , z+ = y/scos { ^^''^ ^ log(g)^ , 

z_ = yssm I log(,s) 1 and W = . 



We define: 



Z±=\ 1 



The assumption <j > yields the bound 

\Z±{s)\ <s for s > 1 (5.75) 

and for all 2 < s < cr: 

|(Z±)i(s)||(Z^)i(c7)| <aloga, |(Z±)i(s)||(Z^)2(a)| <log^, 

|(Z±)2(.)||(Z^)i(a)| < ^loga, \{Z^UsmZ^)2ia)\ < i^. ^^'^^^ 

Moreover, if Z/j solves the homogeneous equation (|5.74p with 

lim Zh = 0, then Zh = 0. (5.77) 



We now turn to the inhomogeneous ordinary differential equation (|5.7ip which we 
solve using the variation of constants method: 

Z{s) = a+{s)Z+{s) + a_(s)Z_(s) (5.78) 

where 

f+^ F,{Z^h-F2{Z^)i 

/+~F2(Z+)i-Fi(Z+)2 ^^■'^^^ 



a. = a_- da, 

where a^. are real constants and where W is the wronskian of z±. Observe from 
()5.72p . (|5.75p . ()5.76p that the integrals in ()5.79p are absolutely convergent, and there 
holds the bounds: 

^7 W ^ Fi{Z-)2-F2{Z^)i . f^°^ F2{Z+)i-F,{Z+)2 
(Z+)i(s) J — dcj + (Z_)i(s) J — da 

r+00 

< / {\Fi{a)\ + \F2{a)\a)logada ^0 as s ^ +00, (5.80) 

J s 
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and similarly: 



< 

s 



Z+h{s) J ^ dcj + (Z_)2(s) J — da 

- (|Fi(cj)| + |F2(cj)|cj)logo-cio- ^ as s ^ +00. (5.81) 

S Js 



We conclude from (l5T7)l that = a° = and (f5T8]l . (I5T9D . (IMD . (ISlSTD now 
yield (|5.73p . This concludes the proof of Lemma 15.41 



Appendix C: Proof of the estimate (15. 4p 

This appendix is devoted to the proof of the estimate (|5.4|) . We in fact claim the 
stronger estimates: 

|A-A,| + |6-6e|<|i^ |a-ae| + |/3-/3c|<|^|^ It - 7c| < (5.82) 
step 1 Improved bound on the modulation equations. 

Let Mod{t) be the vector of modulation equations given by (|3.19p . and rewrite 
the modulation equations (j5.52p using the bounds (j5.3p : 

Mod{t) = n{V) + 0{V'' + V^\\e\\L2) 

= TZ{V) + 0{\tf) (5.83) 

with R(V) polynomial in V satisfying (j5.5ip . Unfortunately, this estimate is just 
not enough to conclude and we need to slightly improve the inner products (|5.53p 
and (j5.55p to gain a cancellation on the null space for the (6, (3) laws from (|5.83p to 

(6, + 62_do(a,a),/3, +6/3-co(q)A-/33A3) = n{V) + O {V^ + V^\\e\\L2) 

= niV) + 0{\t\'^) (5.84) 

with TZ satisfying (j5.5ip . 

Let us assume the improved bound (I5.84p and conclude the proof of (|5.82p . 

We rewrite (j5.83p for M od{t) , M odc{t) using the fact that TZ vanishes at least at 

order 2 at the origin from (I5.5ip : 

\Mod{t) - Modc{t)\ < \t\^\V - Vc\ + 0{\tf). (5.85) 

Using similarly the improved bound (|5.84p . the definition of do (|3.18p and the de- 
generacy a,ac = 0{t^), we get: 

1(6, + 62) -((6,), + 62)1 

+ |(/3, + 6/3 - co(a)A - t^^X") - ((/3e)s + 6^/3^ - co(ac)Ac - /SsAj?)] 
< \t\''\V -Vc\+0{\t{'). (5.86) 
step 2 Estimates for A — Ac and b — be- 

Let us define: 

V = V-Vc = 0{\t\'^) (5.87) 
from (filGOll . (jMH), dlSS]). We claim the bound: 

|6 - bc\ + |A - Acl < \tf + 1^ \V{r)\dT + \t\ (^J^ IH^)!^) • (5.88) 
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Indeed, we first have: 



1 



(5.89) 



Observe now from (|4.6ip . (I4.62p that 

t 

"Co 



A(t) = -^ + 0(|f), b{t) = -^ + 0{\tf) 



from which: 



Injecting this into (|5.89p together with ()5.85p for A and the improved bound (|5.86p 
for b yields: 

Xt-^b + ^X = F^, 6^-^6+Aa = F2, (5.90) 
where Fi and F2 satisfy: 

\Fi\ + \F2\<\tf + \V\. (5.91) 

We rewrite (f5^ 

Zt = -MZ + F, (5.92) 



where Z, M and F are given by: 

A 



Z 



, M 



-1 Co 

Co ^ 



F 



Fx 
F2 



(5.93) 



The eigenvalues of the matrix M are and 1, and hence the system can be rewritten 
in an eigenbasis: 



Zt = —Z + F with D 
t 




1 



(5.94) 



with from HOT]) . (OSD . (ISMIl : 



Z 



\F\<\tt + \Vl j^O as t^O. 



The explicit integration of (|5.94p yields: 

1^(01 < \z{t)\ 



1: 


F{r) 


dr < \t\^+ [ 




r 


Jo 



■\nr)\dT 



\F{T)\dT+\t\ 

which implies ()5.88p . 

step 3 Estimates for a — ac and (3 — fic- 

We claim the bound: 

|a-a,| + |/3-/3e| < |t6ln(|t|)| + |t| 
Indeed, we first have: 

A + ^ _ £oWA±M! = (A + 6^ - .„(„)A - A,A') 
We now derive from (j4.60p - (|4.62p the bounds: 







|log(T) 

\A 



(5.95) 



(5.96) 
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H = + — co(a) + 0(1^1). 



ai + ^^ = Gi, ^,--^+^co(q) = G2, (5.97) 



A2 A2 A2 ' A2 t- t 

Injecting this together with (|5.85|) and the improved bound (|5.86p into (|5.96p yields: 

with: 

\Gi\ + \G2\<t'' + \V\. (5.98) 

We set: 

!? = |, (5.99) 

and rewrite ()5.97p as: 

ai + 2Cor? = Gi, + ^co(a) = ^. 

In order to solve this system, we first diagonalize the matrix cq which is negative 
definite from (j2.13p and hence diagonalizable with eigenvalues ri < 0, r2 < 0. We 
obtain in an eigenbasis: for j = 1,2, 

f (a,)t = -2Cor/. + 0(t5 + P), 



We then perform the change of variables s = and rewrite (jS.lOOp : for j = 1,2 



V 



-2 







\ 2Co— i / \ ^2 

with <f = — TjCq > and 

We have from (jOTj) . (gSZ]), ([S^ST]) the bounds: 

lim Z,- = and F, (s) = O ( 4 ) , 
and we thus conclude from Lemma 15.41 

hjl + < + ^) logad^ < |t'l0g|t|| + p-^^^\V{T)\dT 

which from (j5.99p implies (j5.95p . 

step 4 Bound on P. 

We conclude from l^5M\i . (fOS]) that 

121 < \t'log\t\\ + I' \V{T)\dT + |t| I' \M[A\V{r)\dT. 
Using the initial a priori bound (I5.87P and iterating several times yields: 

\E\<t'\iog\t\\\ 

Injecting this into (|5.88p . (|5.95p yields in return: 

\X-Xc\ + \b-b,\<t^ \a-ac\ + \P-(3c\<\t\^ (5.101) 
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step 5 Bound on the phase parameter. 

We are now in position to derive the key estimate on the phase parameter: 

l7-7c|<t'. (5.102) 

Indeed, we have: 

7t = ^ (7s + («,«)), (5.103) 

and from the improved bound (|5.101|l F^: 

l + |/3|2 + di(a,a) l + \/3c\'^ + di{ac,ac) ^/^l^-^cl,,. .,,1 
1^2 ^2 ^ \ — [ip \- \p - Pc\ + \a - Oc 

= 0{\tf). (5.104) 

dSSSD, (I5J03]) and (ISTM) now yield: 

i, = o{\tf + \v\), 

which after integration in time yields (j5.102p . This concludes the proof of (|5.82p 
assuming (|5.84p . 

step 6 Proof of the improved bound ()5.84p . 

We now turn to the proof of ()5.84p . We change the orthogonality conditions ()3.4p 
and ()3.6p which respectively govern the law of (3 and b for 

(e2,VS + AA^) - (ei,VG) = (5.105) 

and 

(e2,AS + BA2) - (ei,Ae) = (5.106) 
where A and B are well localized functions to be chosen. We claim that for a suitable 
choice A,B, the computation of the modulation equations like for ()5.64p . (|5.69p 
directly leads to the improved bound (|5.84[ ^. 

Indeed, let Mi,M2 be given by (|3.13p . (|3.14p . and let the complex operator for 

e = ei+ie2 ; 

M(e) = Mi(ei,e2) + M2(ei,e2) = Mi(e2, £2) + iM2(ei, £2) - «^Ae + 2i/3 • Ve. 
An integration by parts yields the adjunction formula: 



3fi(M(£),/) = K(£,M(/)). 
Given an orthogonality condition 

9(eJ) = 0, 

the linear term in the computation of the modulation equations is from ()3.1ip (|3.12p 
up to terms 0(t^||£||j;^2): 



-(Mi(£i,£2),/i) + (£2,5,/i) - (M2(£i,£2),/2) " (ei,^,^) = -di{e , M (f) - iOsf) . 
Following closely the proof of ()5.64p . ()5.69p . the improvement ()5.84p follows frorrF^: 

Re{e, M{VQv + AM) - ids{VQv + A2^)) = 0{\tf\\e\\L2), (5.107) 



^^which is critical for the scaling parameter 

We could have chosen this improved set of orthogonality conditions from the very beginning. 

We chose not to do so because this refinement is at first hand not natural and is really needed only 

here, see Remark 15.51 
28,. 



using also the L conservation law H5.63[l 
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Re{e, M{KQr + X^B) - ids{AQv + A^^)) = -2R{e,Q^) + 0{\tf\\e\\L2). (5.108) 
Proof of (|5.1U7p : We rewrite the equation I^THh of Q-p- 

-ib^dbQv + ^Qv -Qv+ ^(^(^)^ + "(^)) q^|q^|2 + i^AQp - 2i/3 • VQ^ 

k[a[t)) 

= ©(t^e-^ls^l). (5.109) 
By differentiating ()5.109p . we obtain: 

-ib^dbVQ-p - Mi(VQp) - iM2{VQv) + ibAVQp - 2i(3 • V^Qp 
= -^^^^^\Qr\'Qr-^bVQr + 0{t^e-\y\) 

which yields: 

MiVQv + X'^A)-idsiVQv + X'^A) = {V^k{0){y,-)Q^ + L+{A)) + ibVQv 

+ Odtpe-I^l). (5.110) 
From (|2.10p . we then chose A solution to 

L,(A) = -V=.(0)(., + ayQ with = 

(VQ,yQ) 

and thus from (|5.110p with / = VQp + A^^: 

Re{e,M{f)-idsf) = b'^ie^VQ^) + aX^^{e,'^) + 0{\tf\\e\\L2) 
= 0{\tf\\e\\L2) 

where we used (|3.5p . ()5.105p in the last step. This concludes the proof of ()5.107p . 
Proof of ()5.108p : We similarly derive from (|5.109p the relation: 

-ib^dbAQv - Mi(AQp) - iM2{AQv) + ibA^Qv - 2i(3 ■ VAQ-p 

= 2ib^dbQp + 2Qv - ^^•^MAj/ + a) |^^|2Q^ _ 2ibAQj, + 2z/3 • VQv + ©(t^g-l?^! 

K{a) 

which implies: 

MiAQv + X'^B) - idsiAQp + A^S) = A^ {V^k{0){y, y)Q^ + L+{B)) 
- 2ib^dbQp-2Q'p + 2ibAQp-2i^-VQ'p + 0{\tfe-\y^). (5.111) 
We hence choose from (|2.10p B solution to: 

L+(i?) = -V2/c(0)(y,y)Q3. 

We now observe that 



dbQv = -i^-^Q + 0{\t\e-\y^] 



so that (|5.11ip yields: 



Re{s, M{AQr + X^B) - ids{AQr + X^B)) + 23f?(e, Qp) 

< h'\R(e,WQ)\ + &|3(e, AQ^)I + /3|9(e, VQ^)| + Odtl'llelU^; 
II^) 



Od^Pllel 



where we used the orthogonality conditions (|5.105p . (|3.7p . (|5.106p . This concludes 

the proof of (fSTTO . 

This concludes the proof of (|5.84p . 
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Remark 5.5. The estimate (j5.82p is stronger than (|5.4p . We need these estimates 
because we use |A — Ad < in the proof of (|5.104p . and thus it is only to control the 
phase paramater in (|5.82p that we need the improved bound (|5.84p . or equivalently 
the refinement of the orthogonality conditions (|3.4p . ()3.6p . 



Appendix D: Algebraic derivation of f lCTD . flOHj) . flOTj) . flOSj) 
Let us recah the definitions l^m . I^m . (lOOll of MW,Mf ^M^ \ 

Proof of (lOi)) . dOH]) : We compute from (lOHl) . (fOOjl . (ICT))) and the control 
of the modulation parameters on Uc- 

M^^\VQ + \lA2)-ids(yQ + XlA2) = L+{VQ + XIA2) - dX^QT^VQ 

- hlV^kiO){y,y)Q^VQ + OiXle-'^\y\) 

where from (|2.12p and the degeneracy \ac\ ^ it suffices to let be the radial 
solution to: 

L+{T^) = ^V^m{y,y)Q\ (5.112) 
To fulfill (|5.34p . we recall L^(VQ) = and chose A2 solution to 

L+{A2) = 6QT^VQ + ^V^k{0){y, y)Q^VQ - awQ (5.113) 

with from (l2Tn]l and the fact that Q, are radial: 

(GQT^d.Q + IV^k{0){y,y)Q^djQ,d,Q) 
«i = Jq2 • (5.114) 

We now observe that ai can be computed in an explicit form. Indeed, we first 
differentiate twice the equation of Q (jl.Sp to derive: 

L+(AQ) = 6|VQpQ, (5.115) 

and hence using the fact that Q, are radial and (j5.112p : for j = 1, 2, 

2{6QT^djQ,djQ) = 6(rO,Q|VQ|2) = (rO,L+(AQ)) = (L+(rO),AQ) 

^v2fc(0)(y,y)g^AQ 

lv^kiO){y,y)Q^VQ,VQ^ - (^^ViV^kmy,y))Q\VQ 

3 
2 



V2fc(0)(y,y)Q2VQ,Vg) +^ j A{V^kmy,y))Q' 
lv'kmy,y)Q'VQ,VQ^ +'-^^^^^ J 
3{V'k{0){y,y)Q'd,Q,d,Q) + '-^^^^ j 



which together with (|5.114p yields: 



tr(V2A;(0)) f , , 

a, = -^^f^ > (5.116) 



since V A;(0) is negative definite, and (|5.36p is proved. 
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Proof of (f05l) : We compute from ([OH]) . dOOjl . (fOOl) and the control of the 
modulation parameters on Uc- 

M^^\i{yQ + XIB2)) - ids [KyQ + A^^a)] = iL^{yQ + A^^a) - 2i\lQT^yQ 

- '-\lV^m{y,y)Q^yQ + 0{Xle-'\y\). 
From L — (yQ) = — 2VQ, (|5.35p follows by choosing B2 solution to: 

= -2A2 + 2Qr°yQ + ^^^^^Q'yQ. 

which is solvable from ()2.1ip since its RHS is orthogonal to Q by radial symmetry 
and dSHSD. 

Proof of (I07)l : 

We compute from ([OH]) . (fOO]) . (fOO]) . (fTTll and the control of the modulation 
parameters on Uc- 

m(^) (AP-p^ + ^2 + ^3 + i^s) - [APv^ +E2 + E3 + iDs] 
+ (co(ac)Ac + ■ y [AP-p^ + E2 + E3 + iD^] 
= M^^\APpJ - idsAT2 + (co(ae)Ae + (i^Xl) ■ yAQ 

+L+{E2) + L+{E3) + iL_(Z)3) - ids{E2) + ©(A^e"!?'!). (5.117) 

We now evaluate M(^)(AP-pJ. From ([231), (I¥:6(1 - (H:62] 1 and the fact that dxPv = 
0{r), daPv = 0(V), dbPp = 0(V^) and a^^P-p = 0{V^), Pv satisfies the following 
equation: 

-iXhdy^Pr + APp - + M^|)^±^p^|p^|2 _ (Aco(a) + ^3A3) . yPp 

= 0(AV'=I^'I). (5.118) 
Differentiating this relation yields: 

ibXdx{AP-p) + M(APp) + (Aco(a) + /SgA^) • 2/APp = -2Pp (5.119) 

+ Ay • ^^(^^ + «) |p^|2p^ _ 3(Aco(a) + /JsA^) • yPp - 2z6A9a(Pp) + ©(A^e'l^l). 

Finally, (|5.117p . (|5.119p and the control of the modulation parameters on Uc yield: 
m(^) [APp^ + E2 + E-i + iD^) - ids [APv, +E2 + E3 + iDs] 
+ (co(ac)Ac + /33A;?) • y [APp^ + + ^3 + ^^3] 

= -idsAT2 - ibXdxAT2 - 2Pp^ + Xy • Y]<^^±Aq^ 

-3(co(ac)A, + (33XI) ■ yQ - 2ibcXcdx{T2) 

+L+{E2) + P+(P3) + iL^{D3) - ids{E2) + O(A^e-l^l) 

= -2Pp, + V^m{y, y)XlQ^ + V^m{y, a,)X,Q^ + ^V^m{y, y, y)XlQ^ 
-3(co(ac)Ac + P3XI) ■ yQ - 2ih,Xcdx{T2) 

+L+{E2) + P+(P3) + iL-{D3) - ids{E2) + 0{Xte-\y\). (5.120) 
To fulfill (|5.37p . we choose E3 and D3 solutions to: 

L+{E2) = -V^k{0){y, y)XlQ^ - V^k{0){y, a,)XcQ^ + 3co{ac)X, ■ yQ + A,a2(ac) • vQ, 

(5.121) 
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where 02 is given by: 



2{V^ k{0){y,.)Q\ djQ) 

iyQ,djQ) 

1. 



(a2(-))i = . ^ o — 3(co(.))i,J = 1,2. (5.122) 



L+(^3) = -^V3A:(0)(y, y, y)\lQ^ + S/SgAj^ • yQ + A^ag • yQ, (5.123) 
where 03 is given by: 

{V^kmy,y,y)Q\d,Q) 

2(yg,9,Q) ^^-^2^) 

and 



L_(I)3) = ^>?cdxT2 - )r\ldx{E2). (5.125) 

Prom ([2^0]) and the definition (|5.122p and (|5.124p of 02 and 03, we may solve (|5.12ip 
and (|5.123p . To conclude the proof of (j5.37p , it remains to prove that we may solve 
(|5.125p . In fact, we have: 

{E2,Q) = -^{L+{E2),AQ) = \{V^m{y,y)Q\KQ)\l = 0, (5.126) 

where we have used (|5.12ip and (|2.16p . Moreover, we have from (j2.27p : 

(r2,Q) = 0. (5.127) 
We deduce from ([2T1|) . (|5.126p and (I5.127P that we may indeed solve (|5.125p . 

Proof of (fOHD : We compute from (fOHjl . (fOOjl . (fOOD and the control of the 
modulation parameters on Uc'- 

M(4)(i(|y|2pp^ +F2)) [i(|y|2pp^ +F2)] 

= M(4)(i|2/|2ppJ + iL_(F2) + 0(A^e-l^l). (5.128) 

We now evaluate M(^)(i|ypP-pJ. We deduce from (|5.118p : 

M(^)(i|2/|2ppJ = -4iAPp^ + 0{\le-\y\). (5.129) 

(I5T:28]1 and (l5J29]l yield: 

M(4)(i(|y|2pp^ + P2)) -ia. [i(|y|2pp, +P2)] 
= -UKPT,^ + iL_{F2) + 0{\le-\y\). 
To fulfill (|5.38p . we chose P2 solution to: 

L_(P2) = -^E2, 
which is possible from (|2.1ip and ()5.126p . 

Proof of (IOQD : (IOQD follows from the fact that M^'^\p) = L+{p) + ©(A^e"'^!) 
and L^{p) = \y\'^Q- 
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